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Abstract 



These lecture notes illustrate the application of Dyson-Schwinger equations in QCD. The exten- 
sive body of work at zero temperature and chemical potential is represented by a selection of 
contemporary studies that focus on solving the Bethe-Salpeter equation, deriving an exact mass 
formula in QCD that describes light and heavy pseudoscalar mesons simultaneously, and the 
calculation of the electromagnetic pion form factor and the vector meson electroproduction cross 
sections. These applications emphasise the qualitative importance of the momentum-dependent 
dressing of elementary Schwinger functions in QCD, which provides a unifying connection be- 
tween disparate phenomena. They provide a solid foundation for an extension of the approach 
to nonzero temperature and chemical potential. The essential, formal elements of this appli- 
cation are described and four contemporary studies employed to exemplify the method and its 
efficacy. They study the demarcation of the phase boundary for deconfinement and chiral sym- 
metry restoration, the calculation of bulk thermodynamic properties of the quark-gluon plasma 
and the response of tt- and /?-meson observables to T and fi. Along the way a continuum order 
parameter for deconfinement is introduced, an anticorrelation between the response of masses 
and decay constants to T and their response to /i elucidated, and a (T, ^)-mirroring of the slow 
approach of bulk thermodynamic quantities to their ultrarelativistic limit highlighted. These 
effects too are tied to the momentum-dependent dressing of the elementary Schwinger functions. 
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1. Introduction. 

In this article 1 describe the application of Dyson-Schwinger equations (DSEs) to QCD 
at finite temperature, T, and quark chemical potential, /i. It is not a pedagogical intro- 
duction, as this can be found in recent reviews. ^'^ The goal instead is to illustrate how 
contemporary studies at (T = 0, = 0) can be used as a foundation and springboard for 
the apphcation of DSEs at finite T and and to describe some of these applications and 
their results. 

The DSEs provide a nonperturbative, Poincare invariant, continuum approach to solving 
quantum field theories. They are an infinite tower of coupled integral equations, with the 
equation for a particular n-point function involving at least one m > n-point function. A 
tractable problem is only obtained if one truncates the system, and historically this has 
provided an impediment to the application of DSEs: a priori it can be difficult to judge 
whether a particular truncation scheme will yield qualitatively or quantitatively reliable 
results for the quantity sought. As integral equations, the analysis of observables using 
DSEs rapidly becomes a numerical problem and hence a critical evaluation of truncation 
schemes often requires, or is at least simplified, by easy access to high-speed computers.^ 
With such tools now commonplace, this evaluation can be pursued fruitfully. 

^The human and computational resources required are still modest compared with those consumed in 
contemporary numerical simulations of lattice-QCD. 
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The development of efficacious truncation schemes is not a purely numerical task, and 
neither is it always obviously systematic. For some, this last point diminishes the appeal 
of the approach. However, with growing community involvement and interest, the quali- 
tatively robust results and intuitive understanding that the DSEs can provide is becoming 
clear. Indeed, someone familiar with the apphcation of DSEs in the late-70s and early-80s 
might be surprised with the progress that has been made. It is now clear ^'^ that trun- 
cations which preserve the global symmetries of a theory; for example, chiral symmetry 
in QCD, are relatively easy to define and implement and, while it is more difficult to 
preserve local gauge symmetries, much progress has been made with Abelian theories^ 
and more is being learnt about non-Abelian ones. 

The simplest truncation scheme for the DSEs is the weak-coupling expansion. Using this 
systematic procedure it is readily established that the DSEs contain perturbation theory, 
in the sense that for any given theory the weak-coupling expansion of the equations 
generates all the diagrams of perturbation theory. Hence, at the very least, the DSEs 
can be used as a generating tool for perturbation theory, and in this application they 
are an essential element in the proof of the renormalisability of a quantum field theory. 
This feature also places a constraint on other truncation schemes; i.e., the scheme must 
ensure that perturbative results are recovered in that domain on which a weak-coupling 
expansion is known to be valid. 

The most important feature of the DSEs is the antithesis of this weak-coupling expan- 
sion: the DSEs are intrinsically nonperturbative. They can be derived directly from 
the generating functional of a quantum field theory and at no stage in this derivation is 
a DSE represented as a sum of diagrams in perturbation theory. Hence their solution 
contains information that is not present in perturbation theory. They are ideal for the 
study of dynamical chiral symmetry breaking^ (DCSB) and confinement in QCD, and of 
hadronic bound state structure and properties. In this application they provide a means 
of elucidating identifiable signatures of the quark-gluon substructure of hadrons. 

Their intrinsically nonperturbative nature also makes them well suited to studying QCD 
at finite-T and yU, where the characteristics of the phase transition to a quark-gluon plasma 
are a primary subject. The order of the transition, the critical exponents, and the response 
of bound states to changes in these intensive variables: all must be elucidated. The latter 
is of particular importance because there lies the signals that will identify the formation 
of the plasma and hence guide the current and future experimental searches. 

There is a significant overlap between contemporary DSE studies and numerical simula- 
tions of lattice-QCD. Of particular importance is that both admit the simultaneous study 
of DCSB and confinement, the absence of which defines the plasma. The DSEs provide an 
adjunct to lattice simulations. They are a means of checking them, and the simulations 
can provide input into the development and constraint of DSE truncations. A truncation 
that is accurate on the common domain can be used to extrapolate into that domain 
presently inaccessible to lattice-simulations, such as finite chemical potential and the T- 
and //-dependence of hadron properties. 

'''Historically, the DSE for a fermion propagator has found widespread use in the study of dynamical 
symmetry breaking; for example, it is the "gap equation" that describes Cooper-pairing in an ordinary 
superconductor. 
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2. Essential Elements of the DSEs. 



In this section I summarise some of the results upon which much of the successful DSE 
phenomenology is founded. Before doing so it is important to specify that I employ a 
Euchdean metric throughout. For real 4-vectors, a, b: 

4 

a • 6 := 6^(5^^ := ^ aj 6j , (1) 

i=l 

and hence a spacelike vector, Qf^, has > 0. The Dirac matrices satisfy 

7i = 7m, = 2 5^^ (2) 

and 75 := -71727374- 

My point of view is that the Euclidean formulation is primary; i.e., a field theory should 
be defined in Euclidean space, where the propagators and vertices are properly called "n- 
point Schwinger functions". This is the perspective adopted in constructive field theory 
and, at least as a pragmatic artifice, by practitioners of lattice-QCD. If the field the- 
ory is well-defined, it is completely specified once all its Schwinger functions are known. 
Analytic continuation in the Euclidean-time variable yields the Wightman functions and, 
following appropriate time-ordering, the Minkowski space propagators. Additional details 
and discussion can be found in Sec. 2.3 of Ref. [1]. 

It is important because the analytic structure of nonperturbatively dressed Schwinger 
functions need not be the same as that of their free-particle seeds. Hence, a priori one 
cannot know the analytic properties of the integrand in a DSE and any rotation of the 
integration contours, as in a "Wick rotation", is plagued by uncertainty: there may be 
poles or branch cuts, etc., that cannot be anticipated from the free-particle form of the 
Schwinger functions involved. This is manifest in the fact that the transcription rules: 

Configuration Space Momentum Space 

d^X^ ^ -ij d'^X^ 1. J d^k"" ^ i J d'k^ 

2. p ^ i7^ -d^ 2. ^ ^ -i7^ • 

3.4^ -i-f^ ■ 3. k^q" -k^ ■ 

4. A^B^" -A^ ■ 4. kf.x^' -k^ ■ x^ , 



are valid at every order in perturbation theory; i.e., the correct Minkowski space integral 
for a given diagram in perturbation theory is obtained by applying these transcription 
rules to the Euclidean integral. However, for skeleton diagrams; i.e., those in which each 
line and vertex represents a fully dressed n-point function, this cannot be guaranteed. 

2.1) Gluon Propagator. 

In Landau gauge the two-point, dressed-gluon Schwinger function, or dressed-gluon 
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Figure 1: DSE for the gluon vacuum polarisation and propagator: sohd hne - quark; 
spring - gluon; dotted-line - ghost. The open circles are irreducible vertices. As indicated, 
the quark loop acts to screen the charge, as in QED, while the gluon loop opposes this. 
Anti-screening the charge and enhancing the interaction. 

propagator, has the form 



where n(A;^) is the ghion vacuum polarisation, which contains all the dynamical infor- 
mation about gluon propagation. This propagator satisfies the DSE illustrated in Fig. 1 
(a nonlinear integral equation). As already stated, a weak-coupling expansion of this 
DSE reproduces perturbation theory. Using this one sees directly that in the one-loop 
expression for the running coupling constant: 



the "llAi'c" comes from the charge- antiscreening gluon loop and the "2A^j" from the 
charge-screening fermion loop, which illustrates how the non-Abelian structure of QCD 
is responsible for asymptotic freedom and suggests that confinement is related to the 
importance of gluon self-interactions. 

Studies of the gluon DSE have been reported by many authors ^ with the conclusion that, 
if the ghost-loop is unimportant, then the charge-antiscreening 3-gluon vertex dominates 
and, relative to the free gauge boson propagator, the dressed gluon propagator is signif- 
icantly enhanced in the vicinity of /c^ = 0. The enhancement persists to ~ 1-2 GeV^, 
where a perturbative analysis becomes quantitatively reliable. In the neighbourhood of 
]<? — ^ the enhancement can be represented^ as a regularisation of as a distribution, 
which is illustrated in Fig. 2. As I will elucidate, a dressed-gluon propagator with the 
illustrated enhancement at /c^ ~ generates confinement and DCSB without fine-tuning. 
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Figure 2: Q{k'^)/k'^ from a solution^ of the gluon DSE (dash-dot hne) compared with 
the one-loop perturbative resuh (dashed hne) and a fit (sohd hne) obtained foHowing the 
method of Ref. [7]; i.e., by requiring that the gluon propagator lead, via the quark DSE, 
to a good description of a range of hadron observables. 



2.2) Quark Propagator. 

In a covariant gauge the two-point, dressed-quark Schwinger function, or dressed-quark 
propagator, can be written in a number of equivalent forms 

Sip) := . (5) 

I'j ■p + S(p) 

-i'y ■ pavip^) + CTsip"^) . (6) 



^7 ■ pA{p'^) + Bip^) 



Tj{p) is the dressed-quark self-energy, which satisfies a nonlinear integral equation: the 
quark DSE (depicted in Fig. 3) 

E(p) = {Z^-l)i-i-p + Z,m^^+ j\^D,,{p-q)^-1^S{q)Vl{q,p), (7) 

where r"(g;p) is the renormalised dressed-quark-gluon vertex, mbm is the A-dependent 
current-quark bare mass that appears in the Lagrangian and := d'^q/{2Tr)^ rep- 
resents mnemonically a translationally-invariant regularisation of the integral, with A 
the regularisation mass-scale. The final stage of any calculation is to remove the regu- 
larisation by taking the limit A oo. The quark-gluon- vertex and quark wave func- 
tion renormalisation constants, Zi(/i^, A^) and Z2(/i^, A^), depend on the renormalisation 
point, /i, and the regularisation mass-scale, as does the mass renormalisation constant 
Z^{^^\^') ■.= Z,{ix',A')-'Z,iix',A'). 
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Figure 3: DSE for the dressed-quark self-energy. The kernel of this equation is constructed 
from the dressed-gluon propagator {D - spring) and the dressed-quark-gluon vertex (F - 
open circle) . One of the vertices is bare (labelled by 7) as required to avoid over-counting. 



One can define a quark mass-function: 

- |g (8) 

and, as depicted in Fig. 4, solving the quark DSE using a dressed-gluon propagator that 
behaves as illustrated in Fig. 2 and a dressed-quark-gluon vertex, r^(p, g), that does 
not exhibit particle-like singularities at {p — g)^ = 0,^ yields a quark mass-function that 
mirrors the infrared enhancement of the dressed-gluon propagator. The results in Fig. 4 
were obtained ^ with the current-quark masses: 

= 3.7MeV, m^ = 82MeV, m^ = 0.59GeV, m^ = 2.0GeV, (9) 

at a renormalisation point of ~ 20GeV. Applying the one-loop evolution formula, 
Eq. (39), these masses correspond to: 

mlfa^^ = 5.5 MeV, m] = 130 MeV, = 1.0 GeV, ml = 3.4 GeV (10) 

and although it is obvious from Fig. 4 that the one-loop formula docs not describe correctly 
the momentum evolution of the mass-function down to = 1 GeV^, the values in Eq. (10) 
provide a useful and meaningful comparison with the values quoted conventionally. 

The quark DSE was also solved in the chiral limit, which in QCD is obtained by setting 
the Lagrangian current-quark bare mass to zero. ^ From the figure one observes immedi- 
ately that the mass-function is nonzero even in this case. That is DCSB: a momentum- 
dependent quark mass generated dynamically in the absence of any term in the action 
that breaks chiral symmetry explicitly. This entails a nonzero value for the quark con- 
densate in the chiral limit. The fact that M(p^) 7^ in the chiral limit is independent of 
the details of the dressed-gluon propagator in Fig. 2; they only affect the magnitude of 
M(p2). 

Figure 4 illustrates that for light quarks {u, d and s) there are two distinct domains: 
perturbative and nonperturbative. In the perturbative domain the magnitude of the 



particle-like singularity is one of the form (P^)~", a G (0, 1]. In this case one can write a spectral 
decomposition for the vertex in which the spectral densities are non-negative. This is impossible if a > 1. 
a = 1 is the ideal case of an isolated, (5-function singularity in the spectral densities and hence an 
isolated, free-particle pole, a G (0, 1) corresponds to an accumulation, at the particle pole, of branch 
points associated with multiparticle production. 
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Figure 4: Dressed-quark mass-function obtained in solving the quark DSE using the 
dressed-gluon propagator of Ref . [7] . 



quark mass-function is governed by the exphcit chiral symmetry breaking mass-scale; 
i.e., the current-quark mass. For < 1 GeV^ the mass-function rises sharply. This is 
the nonperturbative domain where the magnitude of M{p'^) is determined by the DCSB 
mechanism; i.e., the enhancement in the dressed-gluon propagator. This emphasises again 
that DCSB is more than just a nonzero value of the quark condensate in the chiral limit! 
The boundary, at ~ 1 GeV^, is that point where the the enhancement in the dressed- 
gluon propagator becomes significant. 

The solution of — A'P{j?) defines a Euclidean constituent-quark mass, For a given 

quark fiavour, the ratio /m^ is a single, quantitative measure of the importance of the 
DCSB mechanism in modifying the quark's propagation characteristics. As illustrated in 
Eq. (11), obtained using the dressed-gluon propagator in Ref. [7], 
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(11) 



this ratio provides for a natural classification of quarks as either light or heavy. For 
light-quarks the ratio is characteristically 10-100 while for heavy quarks it is only 1-2.^ 
The values of this ratio signal the existence of a characteristic mass-scale associated with 
DCSB, which I will denote by M^. For > the propagation characteristics of a flavour 
with < are altered significantly by the DCSB mechanism, while for flavours with 
^ it is irrelevant, and explicit chiral symmetry breaking dominates. It is apparent 
from the figure that ~ 0.2 GeV~ Aqcd- 

'I'ln my Euclidean metric a true quark mass-pole exhibits itself as a real-p^ solution of + R'P{p^) = 0. 
This is absent in the solutions of the quark DSE illustrated in Fig. 4, which is a manifestation of confine- 
ment, as discussed in Sec. 2.3). 



8 



Figure 5: Illustration of the difFractive electroproduction of a vector meson: e~iV — > e~'Np 
with the transition from photon to vector meson proceeding via a quark loop. The shaded 
region represents a Pomeron-exchange mechanism. 

The effect that the enhancement of the dressed-gluon propagator has on the light-quark 
mass-function is fundamental in QCD and can be identified as the source of many observ- 
able phenomena. Further, that this enhancement little afi^ects heavy-quark propagation 
characteristics at spacelike-p^ provides for many simplifications in the study of heavy- 
meson observables. ^ 

2.3) Confinement. 

One aspect of confinement is the absence of quark and gluon production thresholds in 
colour-singlet-to-singlet iS-matrix amplitudes. This is manifest if, for example, the quark- 
loop illustrated in Fig. 5, which describes ^'^ the difFractive, Pomeron-induced 7 — p 
transition, does not have pinch singularities associated with poles at real-p^ in the quark 
propagators. This is ensured if the dressed-quark and -gluon propagators do not have a 
Lehmann representation. 

What is a Lehmann representation? 

Consider the 2-point free-scalar Schwinger function: 

One can write 

A(.)=rd.4i, (13, 

Jo Z-\- (7 

where in this case the spectral density is 

pix) - — : hm [A(— x — ie) — A(— x -|- ie)] — Sirn^ — x) 

which is non-negative. This is a Lehmann representation: each scalar function necessary 
to completely specify the Schwinger function has a spectral decomposition in which the 
spectral densities are non-negative. Only those functions whose poles or branch points lie 
at timelike, real-A;^ have a Lehmann representation. 

The existence of a Lehmann representation for a dressed-particle propagator is necessary 
if the construction of asymptotic "in" and "out" states for the associated quanta is to 
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Figure 6: Homogeneous Bethe-Salpeter equation for a quark- ant iquark bound state: F is 
the solution, the Bethe-Salpeter amplitude, S is the dressed-quark propagator and K is 
the dressed-quark-antiquark scattering kernel. 



proceed; i.e., it is necessary if these quanta are to propagate to a "detector" . In its absence 
there are no asymptotic states with the quantum numbers of the field whose propagation 
characteristics are described by the Schwinger function. Structurally, the nonexistence 
of a Lehmann representation for the dressed-propagators of elementary fields ensures the 
absence of pinch singularities in loops, such as that illustrated in Fig. 5, and hence the 
absence of quark and gluon production thresholds. 

This mechanism can be generalised and applied to coloured bound states, such as colour- 
antitriplct quark-quark composites (diquarks). In this case a study ^ of the 4-point quark- 
quark scattering matrix shows that it does not have a spectral decomposition with non- 
negative spectral densities and hence there arc no diquark bound states. The same ar- 
gument that demonstrates this absence of diquarks in the spectrum of SU{Nc = 3) also 
proves'' that in SU{Nc = 2) the "baryons", which are necessarily diquarks in this theory, 
are degenerate with the mesons. 

The infrared-enhanced dressed-gluon propagators illustrated in Fig. 2 do not have a 
Lehmann representation. Using forms like this in the kernel of the quark DSE yields 
automatically a dressed-quark 2-point function that does not have a Lehmann represen- 
tation. In this sense confinement breeds confinement, without fine-tuning. 

2.4 ) Hadrons: Bound States. 

In QCD the observed hadrons are composites of the elementary quanta: mesons of a 
quark and antiquark, and baryons of three quarks. Their masses, electromagnetic charge 
radii and other properties can be understood in terms of their substructure by studying 
covariant bound state equations: the Bethe-Salpeter equation (BSE) for mesons and the 
covariant Fadde'ev equation for baryons. 

As a two body problem, the mesons have been studied most extensively. Their internal 
structure is described by a Bcthc-Salpeter amplitude, which is obtained as a solution of 
the homogeneous BSE, illustrated in Fig. 6. For a pseudoscalar bound state the amplitude 
has the form 

iEH{k;P)+-fPFH{k;P) (15) 
+ ^-kk-PGH{k;P) + a^,k^P,HH{k;P)\ , 
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Figure 7: First two sets of contributions to a systematic expansion of the quark- antiquark 
scattering kernel. In this expansion, the propagators are dressed but the vertices are bare. 



where, if the constituents have equal current-quark masses, the scalar functions -E, F, 
G and H are even under k ■ P ^ —k ■ P. In Eq. (15), is a flavour matrix that 
determines the mesonic channel under consideration; e.g., T^^ := (1/2) (A"^ + with 
{X-',j — 1 . . .8} the Gell-Mann matrices. The important new element in the BSE is K, 
the fully-amputated, quark-antiquark scattering kernel: by definition it does not contain 
quark- antiquark to single gauge-boson annihilation diagrams, such as would describe the 
leptonic decay of the pion, nor diagrams that become disconnected by cutting one quark 
and one antiquark line. 

K has a skeleton expansion in terms of the elementary, dressed-particle Schwinger func- 
tions; e.g., the dressed-quark and -gluon propagators. The first two orders in one sys- 
tematic expansion are depicted in Fig. 7. This particular expansion, ^ in concert with 
its analogue for the kernel in the quark DSE, provides a means of constructing a kernel 
that, ordcr-by-ordcr in the number of vertices, ensures the preservation of vector and 
axial-vector Ward-Takahashi identities. This is particularly important in QCD where the 
Goldstone boson nature of the pion must be understood as a consequence of its internal 
structure. 

To proceed with a study of meson properties, one chooses a truncation for K. The homo- 
geneous BSE is then fully specified as a linear integral equation, which is straightforward 
to solve, yielding the bound state mass and amplitude. The "ladder" truncation of K 
combined with the "rainbow" truncation of the quark DSE (F^ 7^ in Fig. 3) is the 
simplest and most often used. The expansion of Fig. 7 allows one to understand why this 
Ward-Takahashi identity preserving truncation is accurate for fiavour-nonsinglet pseu- 
doscalar and vector mesons: there are cancellations between the higher-order diagrams. 
And also why it provides a poor approximation in the study of scalar mesons, where the 
higher-order terms do not cancel, and for fiavour-singlet mesons where it omits timelike 
gluon exchange diagrams. 
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Figure 8: DSE for the axial-vector vertex. The driving term is the bare vertex: i'y^'jfx- 
3. A QCD Mass Formula. 

The chiral-limit axial-vector Ward-Takahashi identity (AV-WTI) 



(16) 



where S :— dia.g{Su, Sd, ■ ■ ■), is the statement of chiral-current conservation in massless 
QCD. It relates the divergence of the inhomogeneous axial- vector vertex, Tf^{k;P), to 
a sum of dressed-quark propagators. The vertex is the solution of the DSE depicted 
in Fig. 8, which involves the quark-antiquark scattering kernel, K. It is therefore not 
surprising that in order to preserve this identity when truncating the DSEs the choice of 
K and the vertex, F^, in the quark DSE, are tied together. The divergence PjuF|^(A;; P) is 
a pseudoscalar and hence contains information about pseudoscalar mesons; i.e., Goldstone 
modes. 



3.1) Dynamical Chiral Symmetry Breaking and Goldstone's Theorem. 
In the chiral limit, the axial- vector vertex has the form^ 



-75 



7^Fr(A;; P) + 7 • kk^Gnik- P) - a^, K HR{k; P) 



(17) 



P., 



+ ff^(A;;P)+ /^-^F^(^;P), 



where: Fr, Or, Hr and are regular as P^ ^ 0; P^t^f,{k; P) ~ 0(p2); THik; P) is the 
pseudoscalar Bethe-Salpeter amphtude in Eq. (15); and the residue of the pseudoscalar 
pole in the axial- vector vertex is fn, the leptonic decay constant: 



IhP, = Z2 f itr \(T'')\,^,S{q+)TH{q;P)S{q. 



(18) 



with the trace over colour, Dirac and flavour indices. This expression is exact: the depen- 
dence of Z2 on the renormalisation point, regularisation mass-scale and gauge parameter 
is just that necessary to ensure that the left-hand-side, Jh, is independent of all these 
things. 



It now follows from the chiral-limit AV-WTI that 



f„E„{k-Q) = Bik-"), 

FR{k;0) + 2fHFH{k;0) = Ai^) , 

GR(k;0) + 2fHGH(k;0) = 2A'(k^), 

HR{k-0) + 2fHHH{k-0) = 0, 



(19) 
(20) 
(21) 
(22) 
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where A^k"^) and B{k'^) are the solutions of the quark DSE in the chiral hmit. As em- 
phasised in Sec. 2.2), the appearance of a B{k'^) ^ solution of the quark DSE in the 
chiral limit signals DCSB. It is an intrinsically nonperturbative effect: in perturbation 
theory B{k'^) oc m, the renormalisation-point independent current-quark mass, and hence 
vanishes in the chiral hmit. Equations (17) and (19)-(22) show that when chiral symmetry 
is dynamically broken: 1) the homogeneous, flavour- nonsinglct, pseudoscalar BSE has a 
massless, = 0, solution; 2) the Bethe-Salpeter amplitude for the massless bound state 
has a term proportional to 75 alone, with the momentum-dependence of Enik^Qi) com- 
pletely determined by that of the scalar part of the quark self energy, in addition to terms 
proportional to other pseudoscalar Dirac structures, Fh, Gh and Hh, that are nonzero 
in general; and 3) the axial- vector vertex, r^(/c;P), is dominated by the pseudoscalar 
bound state pole for ~ 0. The converse is also true. 

Hence, in the chiral limit, the pion is a massless composite of a quark and an antiquark, 
each of which has an effective mass ~ 450 MeV. With a dressed-gluon propagator of 
the type depicted in Fig. 2, this occurs without fine-tuning. 



3.2) Nonzero Quark Masses: A Mass Formula. 
When the current-quark masses are nonzero the AV-WTI is modified: 

-iP^T^^ik- P) = 5-^(A;+)75^ + l,^S-\k^) - M(^) rf (A;; P) - T^{k- P) M(^) , (23) 

where: M(^) = diag(m(^, m^, m^, . . .), is the current-quark mass matrix. In this case both 
the axial- vector and the pseudoscalar vertices have a pseudoscalar pole: i.e., 

^U^-.P) = ^l,[l^Fi{k-P)+^-kk,G"^{k-P)-a,,k,H^{k-P)\ (24) 



+ ffr^;P)+/^ 



P'^ + m^H 



rH{k;P), 



and 



r^{k;P) = —j,[iSi{k;P)+^-PJ^i{k;P)+^-kk-Pg^{k;P) 

1 



(25) 



+a^,k^P,ni{k;P) 



P^ + mjj 



TH{k-P), 



with: £f , , Gf , g^, Hg, nl and fg, regular as P^ ^ -m\ and P^V^^{k- P) 

0(P2). The AV-WTI entails^ that 



(m) 



(26) 



where fn is given by Eq. (18), with massive quark propagators in this case, and the 
residue of the pole in the pseudoscalar vertex is 



{T^)\,S{q+)TH{q;P)S{q. 



(27) 



The renormalisation constant Z4 on the right-hand-side depends on the gauge parameter, 
the regularisation mass-scale and the renormalisation point. This dependence is exactly 
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that required to ensure that: 1) is finite in the hmit A — > oo; 2) Th is gauge-parameter 
independent; and 3) the right- hand- side of Eq. (26) is renormahsation point independent. 
This is obvious at one-loop order, especially in Landau-gauge where Z2 = 1 and hence 

= Zjn- 

Equation (26) is a mass formula for flavour-octet pseudoscalar mesons that is valid in- 
dependent of the magnitude of the current-quark masses of a meson's constituents. For 
small current-quark masses, using Eqs. (15) and (19)-(22), Eq. (27) yields 

= - 4r > - - Ul^', A') r troi^ac [Srn=^{q)] , (28) 

Jh •'i 

where the superscript "0" denotes that the quantity is evaluated in the chiral limit and 
, as defined here, is the chiral limit vacuum quark condensate., which is renormalisation- 
point dependent but independent of the gauge parameter and the regularisation mass- 
scale. Hence Eq. (25) is the statement that in the chiral limit the residue of the hound 
state pole in the flavour-nonsinglet pseudoscalar vertex is (— {qq)1)/ fn- 

Now one obtains immediately from Eqs. (26) and (28) 



f>il = -[< + <]{m + 0{ml) (29) 

fl^mU = - K + m^](?-g)J + 0(m9 , (30) 

which exemplify what is commonly known as the Gcll-Mann-Oakcs-Rcnncr relation. 

The primary result, Eq. (26), is valid independent of the magnitude of rhg, and can be 
rewritten in the form 

fnTul^ -{qq)^MH (31) 
where I have introduced the notation 



- (qq)^ fnTH = fnZ, f \ir ^{t")' ^^S{q+)TH{q;P)S{q_) 



(32) 



which defines an in-meson condensate. This emphasises that, for nonzero current-quark 
masses, Eq. (26) does not involve a difference of vacuum massive-quark condensates; a 
phenomenological assumption often employed. 

As elucidated elsewhere, Eq. (26) has another important corollary: it predicts that the 
mass of a heavy pseudoscalar meson rises linearly with the current-quark mass of its heavy 
constituent (s). Model calculations^ show that this hnear evolution is dominant at ^ 2ms, 
in agreement with experiment where the mass of the K, D and B mesons lie precisely on 
the same straight line. 

In Eq. (26) one therefore has a single mass formula, exact in QCD, that provides a unified 
understanding of light- and heavy-meson masses. 



4. An Illustrative Model. 

I have already made use of a model ^ in illustrating some of the robust results of DSE 
studies. To further elucidate the method I will describe that model in more detail. For 
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the kernel of the quark DSE it employs the analogue of the lowest-order BSE kernel in 
Fig. 7: 

j\'D,^{p-q)^^,Siq)n{q,p)^ j\{{p ' <lf) D^P ' <l)^l,S{q)^^, . (33) 

This is the "rainbow" approximation, in which the specification of the model is complete 
once a form is chosen for the "effective coupling" Q {k"^) . 



A choice for Qik"^) can be motivated by observing that at large-Q^ : = 
asymptotically free theory the quark-antiquark scattering kernel satisfies 

g\ll')D,,{p-q) [Tl{p+,q+)S{q+)\ x [S{qJ)Tl{q_,p_) 



{p - qf 



m an 



(34) 



47ra(g2)D;:-(p-g) 



X 



A"" 



where P is the total momentum of the quark-antiquark pair, p+ :— p-\- rjpP and p_ := 
p — {l — rip)P with < 77 < 1. Choosing a truncation of K in which this right-hand-side is 
identified with the lowest order contribution in Fig. 7 then consistency with the AV-WTI 
requires 

GiQ^) ■.= AnaiQ^). (35) 

Thus the form of Q{Q^) at large-Q^ is fixed by that of the running coupling constant. 
This Ansatz is often described as the "Abelian approximation" because the left- and 
right-hand-sides are equal in QED. In QCD, equality between the two sides of Eq. (35) 
cannot be obtained easily by a selective resummation of diagrams. As reviewed in Ref. [1], 
Eqs. (5.1) to (5.8), it can only be achieved by enforcing equality between the renormali- 
sation constants for the ghost-gluon vertex and ghost wave function: Zi = Z^. 



The explicit form of the Ansatz employed in Ref. [7] is 



A;2 



87r^D5\k) + ^D^e-'' + Air ■ 



In 



r+(l + kyAlcBy 



QCD 



Hk') , (36) 



0.234 GeV. 



with J^{P) - [1 - exp(-A;V[4m2])]/P and r = e^ - l,Nf = 4 and A 

The qualitative features of Eq. (36) are clear. The first term is an integrable infrared 
singularity^^ and the second is a finite-width approximation to 5'^(A;), normalised such 
that it has the same / d'^k as the first term. In this way the infrared singularity is split 
into the sum of a zero- width and a finite- width piece. The last term in Eq. (36) is 
proportional to a{k^)/k'^ at large spacelike- A;^ and has no singularity on the real-/c^ axis. 
Gluon confinement is manifest since Q (/c^) / k'^ doesn't have a Lehmann representation. 



4-1) Solving the Quark DSE. 

There are ostensibly three parameters in Eq. (36): D, uj and m^. However, in the 
numerical studies the values u — 0.3GeV(= l/[.66fm]) and rrit — 0.5GeV(= l/[.39fm]) 
were fixed, and only D and the renormalised u/d- and s-current-quark masses varied in 
order to satisfy the goal of a good description of low-energy tt- and X-meson properties. 
This was achieved with 

D = 0.781 GeV^ m(;,^ = 3.74MeV, m^ = 82.5MeV , (37) 
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Figure 9: Ansatz for Q{k'^)/k'^ employed in Ref. [7]. "B&P" labels a solution^ of the gluon 
DSE, which is presented for comparison, as is the one-loop running coupling in QCD. 



at /X ~ 20 GeV, which is large enough to be in the perturbative domain. The effective 
coupling obtained is depicted in Fig. 9. 



Using Eqs. (33), (36), (37), and the renormalisation boundary condition 



i'y ■ p + 



(St 



the quark DSE, Eq. (7), is completely specified and can be solved by iteration. 



The chiral limit in QCD is unambiguously defined by m = 0. In this case there is no 
perturbative contribution to the scalar piece of the quark self energy, _B(p^,yU^): in fact, 
there is no scalar, mass-like divergence in the perturbative calculation of the self energy. It 
follows that Z2{fJ,'^, A^) mbm(A^) = , VA and, from Eq. (38), that there is no subtraction in 
the equation for B{p'^, /x^). In terms of the renormalised current-quark mass the existence 
of DCSB means that, in the chiral limit, M(yU^) ~ 0(l//i^), up to In/i^-corrections.^ 

Figure 10 depicts the renormalised dressed-quark mass function, M(p^), obtained by 
solving the quark DSE using the parameters in Eq. (37), and in the chiral limit. It 
is a complement to Fig. 4 because it highlights the qualitative difference between the 
behaviour of M(p^) in the chiral limit and in the presence of explicit chiral symmetry 
breaking. In the latter case 



M(p2 



large— 


m 




(iln 


p2 




^^QCD 





{1 + two loop} , 7„ 
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33 - 2A^ 



(39) 



TThis is a model-independent statement; i.e., it is true in any study that preserves the one-loop 
renormalisation group behaviour of QCD. 
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Figure 10: The renormalised dressed-quark mass function, M(p^), obtained in solving 
the quark DSE using the parameters in Eq. (37): u/d-qnenk (sohd hne); s-quark (long- 
dashed line); and chiral limit (dashed line). The renormalisation point is yU 20 GeV. The 
intersection of the line M^(p^) = (dotted line) with each curve defines the Euclidean 
constituent- quark mass, . 



However, in the chiral limit the ultraviolet behaviour is given by 

(- (mr) 



M(p^) ^"'^^p' ^^ihn 



In 



7? — 

/'■QCD 



l-7rr 



(40) 



where is the renormalisation-point-independent vacuum quark condensate. Analysing 
the chiral limit solution yields 

= (0.227 GeV)^ (41) 

which is a reliable means of determining (qq)^ because corrections to Eq. (40) are sup- 
pressed by powers of K^qj^/ [x^. 

Equation (28) defines the renormalisation-point-dependent vacuum quark condensate 

^=19 GeV I [Sni=o{q)] 



= (0.275 GeV)= 



/i=19 GeV 



(42) 

the calculated result. It is straightforward to establish explicitly that {qq)1i = constant, 
independent of /i, and hence 

m^{qq)l:^m{qq)\ (43) 

which unambiguously defines the renormalisation-point- independent current-quark masses. 
Prom this and Eqs. (37), (41) and (42) one obtains the values of these masses appropriate 
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to this model 

m„/d = 6.60 MeV , rhs = 147 MeV . (44) 

Using the one-loop evolution in Eq. (39) these values yield m'l^^^ = 3.2 MeV and = 
72 MeV, which are within ~ 10% of the actual values in Eq. (37). This indicates that 
higher-loop corrections to the one-loop formulae, which are present in the solution of the 
integral equation as made evident hy A{p^ , fj,'^) ^ 1 in Landau gauge, provide contributions 
of < 10% at — ji^. The higher-loop contributions decrease with increasing p^. 

Prom the renormalisation-point-invariant product in Eq. (43) one obtains 

- (?"?)Jt=iGeV '-^ [I/Aqcd])^- = (0.241 GeVf . (45) 

This result can be compared directly with the value of the quark condensate employed in 
contemporary phenomenological studies: (0.236 ± 0.008 GeV)^. Increasing a; — > 1.5 a; in 
^(/c^) raises the calculated value in Eq. (45) by ~ 10%, a weak sensitivity. 

After this discussion of the vacuum quark condensate it is straightforward to determine 
the accuracy of Eqs. (29) and (30). Using experimental values on the left-hand-side, one 
finds: 

(0.0924 X 0.1385)2 = (0.113 GeV)^ cf . (0.111 GeV)^ = 2 x 0.0055 x 0.24^ (46) 
(0.113 X 0.495)2 = (0.237 GeV)^ cf. (0.206 GeV)^ = (0.0055 + 0.13) x 0.24^ (47) 

which indicates that 0(m2)-corrections begin to become important at current-quark masses 
near that of the s-quark, as demonstrated further in Ref . [8] . 

4-2) Solving the Pseudoscalar Meson BSE. 

The model quark DSE described above employs the rainbow truncation. Following Fig. 7 
the consistent Ward-Takahashi identity preserving truncation of the quark-antiquark scat- 
tering kernel is the ladder approximation: 

KZiQ. k;P)-- g{{k - qf) D^:r{k - q) , (48) 

\ / tr \ / su 

in which case the explicit form of Fig. (6) is 

VH{k- P) + g{{k - q f) D^;::%k - q)^^^S{q^)rH{q; P)S{q.)^^. = . (49) 

Having an Ansatzior ^(fc^), S{q) in Eq. (49) follows by solving the quark DSE. The kernel 
of the BSE is then completely specified and solving the equation for TH{k\P) and the 
bound state mass is a straightforward numerical problem. Then, with Df^,^{k), S{p) and 
Tnik'jP), the calculation of other observables such as: the leptonic decay constant, fn] 
meson charge radius, th', and electromagnetic form factor, Fh{Q'^)', etc., is possible. 

The general form of the solution of Eq. (49) is given in Eq. (15), where the scalar functions 
depend on the variables k'^ and k ■ P and are labelled by the eigenvalue P^. From this 
it is clear that the integrand in Eq. (49) depends on the scalars: A;^, k ■ q, q^, q - P and 
P"^, which takes a fixed-value at the solution; i.e., at each value of kernel is a 
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All amplitudes 


IT 


chiral limit 
mo /° TIh 


ss 

friss fss T^H 


Method (A) 

Uo only 
Uq and Ui 


0.1385 0.0924 1.01 

0.136 0.0999 0.95 
0.1385 0.0925 1.00 


0.0 0.0898 1.00 

0.0 0.0972 0.94 
0.0 0.0898 1.00 


0.685 0.129 1.00 

0.675 0.137 0.95 
0.685 0.129 1.00 


E only 








Method (A) 
Uq only 


0.105 0.0667 1.82 
0.105 0.0667 1.82 


0.0 0.0649 1.81 
0.0 0.0649 1.81 


0.512 0.092 1.68 
0.513 0.092 1.69 


E, F 








Method (A) 
Uq only 


0.136 0.0992 0.95 
0.136 0.0992 0.95 


0.0 0.0965 0.95 
0.0 0.0965 0.95 


0.677 0.137 0.95 
0.678 0.138 0.95 


E, F, G 








Method (A) 
Uq only 
Uo and Ui 


0.140 0.0917 1.01 
0.136 0.0992 0.95 
0.140 0.0917 1.01 


0.0 0.0891 1.00 
0.0 0.0965 0.95 
0.0 0.0891 1.00 


0.688 0.128 1.01 
0.678 0.138 0.95 
0.689 0.128 1.01 



Table 1: Calculated values of the properties of light, pseudoscalar mesons composed of 
a quark and antiquark of equal- mass. The mass {m^ — 0.1385) and decay constant 
{f^^^ = 0.0924) are in GeV, TZh is dimensionless. With the exception of the calculations 
that retain only the zeroth Chcbyshev moment, labelled by "f/o only", the results are 
independent of the momentum partitioning parameter, f]p, in Eq. (34). 

function of four, independent variables. Solving Eq. (49) can therefore require large-scale 
computing resources, especially since there are four, independent scalar functions in the 
general form of the solution. 

Two different solution techniques can be employed. In one procedure, which I will label: 
(A), the scalar functions are treated directly as dependent on two, independent variables: 
E^k"^, k ■ P: P"^). etc. This requires straightforward, multidimensional integration at ev- 
ery iteration. Storing the multidimensional kernel requires a large amount of computer 
memory but the iteration proceeds quickly. 

An adjunct, which I will label: (B), employs a Chebyshev decomposition procedure. It is 
implemented by writing 

ax 

E{k^,k- P-P^)^ J2 'E{k'^;P^)Ui{cos(3), (50) 

with similar expansions for F, G := k ■ P G and i7, where k ■ P := cos (5\fWP^ and 
{Ui{x)]i = 0, . . . , oo} are Chebyshev polynomials of the second kind, orthonormalised 
according to: 

- I dx Vl -x2 Ui{x)UAx) = Sij . (51) 

TT J-1 

This procedure requires a large amount of time to set up the kernel but does not require 
large amounts of computer memory. 
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All A m nl ifn npc; 


71 j-i 
'IP 


— 50 




71 r> 
'IP 


— 25 




71 n 
'IP 


— 00 

— yj . yjyj 






rriK 


Ik 


TZk 


rriK 


Jk 


TZk 


rriK 


fx 


TZk 


Method (A) 


0.497 


0.109 


1.01 


0.497 


0.109 


1.01 


0.497 


0.109 


1.01 


Uq only 


0.469 


0.117 


0.96 


0.482 


0.117 


0.95 


0.475 


0.113 


1.02 


Uo and Ui 


0.500 


0.111 


1.00 


0.497 


0.109 


1.01 


0.498 


0.110 


1.00 


Uq, Ui and U2 


0.497 


0.109 


1.01 


0.497 


0.109 


1.01 


0.496 


0.109 


1.01 


E only 








Method (A) 


0.430 


0.079 


1.55 


0.430 


0.079 


1.55 


0.429 


0.076 


1.55 


Uq only 


0.380 


0.077 


1.54 


0.401 


0.076 


1.51 


0.415 


0.073 


1.55 


Uo and Ui 


0.439 


0.089 


1.52 


0.430 


0.078 


1.55 


0.431 


0.076 


1.57 


Uo, Ui and U2 


0.430 


0.078 


1.55 


0.430 


0.078 


1.55 


0.427 


0.076 


1.55 


E, F 








Method (A) 


0.587 


0.17 


0.79 


0.557 


0.14 


0.86 


0.533 


0.11 


0.94 


Uo only 


0.505 


0.12 


0.82 


0.518 


0.11 


0.86 


0.512 


0.11 


0.96 


Uq and Ui 








0.556 


0.14 


0.86 


0.537 


0.12 


0.94 


Uq, Ui and U2 


0.583 


0.16 


0.79 


0.557 


0.14 


0.86 


0.532 


0.12 


0.93 


E, F, G 








Method (A) 


0.500 


0.108 


1.01 


0.500 


0.108 


1.01 


0.500 


0.108 


1.01 


Uo only 


0.471 


0.116 


0.96 


0.484 


0.116 


0.95 


0.477 


0.112 


1.02 


Uo and Ui 


0.504 


0.110 


1.00 


0.500 


0.108 


1.01 


0.502 


0.109 


1.00 


Uo, Ui and U2 


0.500 


0.108 


1.01 


0.500 


0.108 


1.01 


0.499 


0.108 


1.01 



Table 2: Calculated properties of the X-meson for various values of the momentum par- 
titioning parameter, r]p; "— " means that no bound state solution exists in this case. The 
mass {rrfK^ — 0.496) and decay constant (/^^ = 0.113) are in GeV, TZk is dimensionless. 



In Tables 1 and 2 I list values of the dimensionless ratio 

i^H — ,20 ■ (52) 

Jh^h 

A value of TZh = 1 means that Eq. (26) is satisfied and hence so is the AV-WTI.t Looking 
at the tabulated values of TZh it is clear that the scalar function H is not quantitatively 
important, with the AV-WTI being satisfied numerically with the retention of E, F and 
G in the pseudoscalar meson Bethe-Salpeter amplitude. The values of TZh, and the other 
tabulated quantities, highlight the importance of F and G: F is the most important of 
these functions but G nevertheless provides a significant contribution, particularly for 
bound states of unequal-mass constituents. 



^It illustrates that the pseudoscalar- meson pole in the axial- vector vertex is related to the pseudoscalar- 
meson pole in the pseudoscalar vertex in the manner elucidated above. A finite value in the chiral limit 
emphasises that oc A4h as A4h 0. 
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Figure 11: The zeroth Chebyshev moment of the scalar functions in the mesons' Bethe- 
Salpeter amphtude: chiral hmit (dotted hne) ; 7r-meson (sohd hne) ; X-meson (long-dashed 
line); and fictitious, ss bound state (dashed line). For ease of comparison the amplitudes 
are all rescaled so that ^Enik"^ = 0) = 1. 

From Tables 1 and 2, and Eqs. (31), (32), (37) and (45), one calculates 

-{m)l=lGeV -{m)^=lGey -{m)%lGeV 

(53) 

(0.245 GeV)3 (0.284 GeV)^ (0.317 GeV)=^ 

showing that, for light pseudoscalars, the "in- meson condensate" 1 have defined increases 
with increasing bound state mass; as does the leptonic decay constant, f^-^ Both of these 
trends are modified as one moves into the heavy-quark domain; —{qq}^ const, and 
^ as Mh oo. 

Figure 11 depicts the scalar functions in the Bethe-Salpeter amplitude obtained as solu- 
tions of Eq. (49), focusing on the zeroth Chebyshev moment of each function, which is 

H~{qq)^)/ Ih is the residue of the bound state pole in the pseudoscalar vertex, just as /h is the 
residue of the bound state pole in the axial-vector vertex. As expected, (99)J^=iGeV ~ (QI)^] -iceV" 
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k' (GeV^) 



Figure 12: Asymptotic behaviour of the 0th Chebyshev moments of the functions in the 
TT-meson Bethe-Salpeter ampUtude: /^°E^(/c^) (GeV, soUd hne); /^'^^^(A;^) (dimension- 
less, long-dashed line); A;^ f 7r(/i^^) (dimensionless, dashed line); and /c^ °iJ^(A;^) (GeV, 
dot-dashed line). The /c^-dependence is identical to that of the chiral-limit quark mass 
function, M(p^), Eq. (40). For other pseudoscalar mesons the momentum dependence of 
these functions is qualitatively the same, although the normalising magnitude differs. 

obtained via 

^EHik"^) ■= - r d(3 sin^ (3 C/o(cos (3) Enik'^, k ■ P; P^) , (54) 

vr JO 

and similarly for F, G [G for the /T-meson] and H. I note that: the momentum-space 
width of ^EHik"^) increases as the current-quark mass of the bound state constituents 
increases; ^Fnik'^ — 0) decreases with increasing current-quark mass but that ^Fnik"^) is 
still larger at k"^ > 0.5 GeV^ for bound states of higher mass; ^nik"^) [^K{k'^)] behaves 
similarly; and the same is true for Huik; P) but it is uniformly small in magnitude thereby 
explaining its quantitative insignificance. 

Figure 12 depicts the large-A;^ behaviour of the scalar functions in the pseudoscalar Bethe- 
Salpeter amplitude. The momentum dependence of ^EH{k'^) at large-A;^ is identical to that 
of the chiral-hmit quark mass function, M{p^) in Eq. (40),^"^ and characterises the form 
of the quark-quark interaction in the ultraviolet. Figure 12 elucidates that this is also 
true of °Ff^(fc2), k^Vnik'^) [k'^^Ki.k'^) for the TT-meson] and k'^^Hnik'^). Each of these 
functions reaches its ultraviolet limit by k"^ ^ 10 GeV^, which is very-much- less-than the 
renormalisation point, /^^ = 361 GeV^. As I will illustrate below, this result has important 
implications for the behaviour of pseudoscalar meson form factors. 

A direct verification of Eqs. (19)-(22) is possible in this concrete model. Consider the 
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Figure 13: An illustration of the realisation of the identities Eqs. (19) and (20), which are 
a necessary consequence of preserving the axial vector Ward-Takahashi identity. 

inhomogeneous axial-vector vertex equation, Fig. 8, in the ladder truncation: 

rf,(A:;P) = ^2757.^-^^^((fc-g)')^;r(^-g)y7M5(g+)rg,(g;P)5(g_)y7.. (55) 

From the homogeneous BSE one already has the equations satisfied by Eu{k\ 0), Fuik; 0), 
Gnik^O) and Hnik^O). To proceed, one substitutes Eq. (17) for T^^{k; P) in Eq. (55). 
Using the coupled equations for Enik; 0), etc., one can identify and eliminate each of the 
pole terms associated with the pseudoscalar bound state. [That the homogeneous BSE 
is linear in the Bethe-Salpeter amplitude allows this.] That yields a system of coupled 
equations for Fr(/c; 0), G_r(/c; 0) and Hji{k; 0), which can be solved without complication. 
[The factor of Z2 automatically ensures that F^lk"^ = /x^; P = 0) = 1.] The realisation of 
the first two identities, Eqs. (19) and (20), is illustrated in Fig. 13. The remaining two 
identities, Eqs. (21) and (22), are realised in a similar fashion. 

5. Additional Phenomenological Applications. 

In the model illustration of Sec. 4) an algebraic Ansatz was developed for the dressed- 
gluon propagator based on the qualitative behaviour of solutions of the gluon DSE. From 
this basic beginning, I illustrated how one can proceed to calculate hadronic observables. 
A number of qualitatively significant features emerged, among them DCSB and confine- 
ment, all of which are related to the strong momentum dependence of the quark mass 
function, Eq. (8), in the infrared. 

That last observation suggests an alternative phenomenological approach: develop an 
algebraic Ansatz for the dressed-quark propagator. This is not as fundamental as the 
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approach in Sec. 4) because many, apparently distinct features of the dressed-quark prop- 
agator are encoded in a few parameters of the Ansatz for the dressed-gluon propagator; 
modeUing the dressed-quark propagator directly requires additional parameters to de- 
scribe correlated effects. However, it has a significant merit: with an algebraic as opposed 
to a numerical representation of the dressed-quark propagator it is possible to study scat- 
tering observables more quickly and easily. The approach can yield quantitatively reliable 
results provided the Ansatz exhibits those essential qualitative features manifest in a di- 
rect solution of the quark DSE using a realistic Ansatz for the dressed-gluon propagator. 
Further, in allowing a rapid analysis of a broad range of observable phenomena, it provides 
a means of exploring the hypothesis that the bulk of hadronic phenomena are simply a 
manifestation of the nonperturbative, momentum-dependent dressing of the elementary 
Schwinger functions in QCD. 



A simple and widely used model is 
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as{i) = 2mJ'{2{i + m'))+r{hi):F{hi){ho + h2Hei)) , (56) 

^^^^^ = 2{^ + m-f ' ^^^^ 

with (Tsid •= ^(^siv^)-, ^v(0 '■— -^"^ (^v{p'^): where = A^.^, A is a mass-scale, and 
J-^{y) := [1 — exp{—y)]/y. This five-parameter form, where fh is the current-quark mass, 
combines the effects of confinement''^ and DCSB with free-particle behaviour at large, 
spacehke p^} Applying Eq. (28) in this case: 

3 r^^ 

-{qq)l := Jim^ Z,{^^' , M') — dssa^s), (58) 
Aqcd 47r2 bih 



and the pion mass is given by 



mlf^^2m{qq)'l^,y. (60) 



In Sec. 4.1) we saw that when all the components of are retained, Eq. (60) yields an 
approximation to the pion mass found in a solution of the Bethe-Salpeter equation that 
is accurate to 2%. 

The model has been used for both u/d- and s-quark propagators with the difference be- 
tween fiavours manifest in 6q b^^'^, 6| ^ b2^^ and ^ rriu/d' the first allows for a 
difference between the K and vr in-meson condensates and the second for Mf ^ ^u/d^ 
and all three arc phenomenological constraints observed in the previous section. As em- 
phasised above, in a solution of the quark DSE using an Ansatz for the dressed-gluon 



■I^The representation of S{p) as an entire function is motivated by the algebraic solutions of Eq. (7) in 
Refs. [16,17]. 

^At largc-j;^: (Tv (p^) ^ l/p^ and crs(p^) ~ m/p^. Therefore the parametrisation does not incorporate 
the additional Inp^-suppression characteristic of QCD: it corresponds to 7m ^ 1 in Eq. (40). It is a 
useful but not necessary simplification, which introduces model artefacts that are easily identified and 
accounted for. e = 0.01 is introduced only to decouple the large- and intermediate-p^ domains. 
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Impulse Approximation 




P 

Figure 14: Impulse approximation to F.„{q^): S labels the dressed-quark propagator; 
the pion Bethe-Salpeter amplitudes; and the dressed-quark-photon vertex. 



propagator, the parameters in Eq. (56) are correlated and one can anticipate this crudely 
when fitting them. 

5.1) Pion Electromagnetic Form Factor. 
The renormalised impulse approximation to the electromagnetic pion form factor is 



(Pi +P2),xF^{q^) := A^(pi,P2) 



(61) 



2Nr 



I 



tVD 



Gnik; -p2)S{k++) iTl{k++, k+-)S{k+-) Q^{k - q/2;pi)S{k. 



kaf3 ■— k + api/2 + /3q/2 and p2 :— pi + q, illustrated in Fig. 14. Here QT^{k] P) is the 
pion Bcthc-Salpeter amplitude normalised such that i?(0; 0) = -B(O) in the chiral limit, in 
which case the consistent canonical normalisation of the Bethe-Salpeter amplitude is 



25'^NlP^ 



tr 



+tr 



dP,.. 



Ql{q--P)S{q+)gi{q-P) 



(62) 



dS{q^) 
dp. 



where QiriQ^—PY '■= C ^ Q.„{—q; —P) C with C = 7274, the charge conjugation matrix, 
and X* denotes the matrix transpose of X. 

Given S and Eqs. (19)-(22), what form does the Bcthc-Salpctcr amplitude take? 

In Sec. 4.2) wc saw that the zcroth Chcbyshcv moments of the pion Bcthc-Salpctcr am- 
plitude provided results for and f.^ that were indistinguishable from those obtained 
with the full solution. Also H^^ ~ and hence it was quantitatively unimportant in the 
calculation of m,r and /^r. These results are not specific to that particular model; in the 
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latter case because the right-hand-side of Eq. (22) is zero and hence in general there is no 
"seed" term for Ht,. We also saw that at large- /c^, independent of assumptions about the 
form of K, one has 

E°(^2;Ooc (63) 

and that the same is true of F^{k'^]P'^), k"^ Gl(k'^; P^) and k'^ H^{k'^; P'^). This makes 
manifest the "hard-gluon" contribution to F^^lq^) in Eq. (61). In addition, in an asymp- 
totically free theory, one has ^ 

k' Gl(k'; P') = 2F^(k'; P') , k' > M^^ , (64) 

with Muv := 10 Aqcd- 

These observations, combined with Eqs. (19)-(22), motivate a model for Qt^: 

E4k;P) = Bo{e) (65) 

with F^{k;P) = E^{k; P)/ {110 f^), G^{k;P) = 2F^{k;p)/[k^ + M^^] and H^ik; P) = 0. 
The relative magnitude of these functions at large k'^ is chosen to reproduce the numerical 
results in Fig. 12. 

The final element in Eq. (61) is r^(pi,p2), the renormalised, dressed-quark-photon vertex, 
and it is because this vertex must satisfy the vector Ward-Takahashi identity: 

{pi-P2),tri{pi,P2) = S-\p^) - S-\p2) , (66) 

that (pi — P2)/i A^(Pi,P2) = and no renormalisation constants appear explicitly in 
Eq. (61). TJ^{pi,p2) has been much studied^ and, although its exact form remains un- 
known, its robust qualitative features have been elucidated so that a phenomenologically 
efficacious Ansatz has emerged 



i^Z{P,Q) — ^^Aip'^,q'^)7^^+ {p + q)^ 



-i7 -{p + q) ^a{p\ r) + ^b{p\ r) 



;(67) 



S/(/,g^) := i[/(p^) + /(g^)], A,(p^g2):=:^^^l-^, (68) 

where f = A, B. A feature of Eq. (67) is that the vertex is completely determined by 
the renormalised dressed-quark propagator. In Landau gauge the quantitative effect of 
modifications, such as that canvassed in Ref. [19], is small and can be compensated for 
by small changes in the parameters that characterise a given model study. 

The model parameters were determined by optimising a least-squares fit to and 
(??)iGeV' ^^<i ^ selection of pion form factor data on the domain q'^ e [0,4] GeV^. The 
procedure does not yield a unique parameter set with, for example, the two sets: 

A(GeV) fh bo h &2 h 
A 0.473 0.0127 0.329 1.51 0.429 0.430, (69) 
B 0.484 0.0125 0.314 1.63 0.445 0.405, 

providing equally good fits, as illustrated in Table 3. There is a domain of parameter 
sets that satisfy the fitting criterion and they are distinguished only by the calculated 
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Calculated Experiment 



0.092 GeV 0.092 

(-(??)iGev)^^^ 0.236 0.236 ± 0.008^3 

ruu/d 0.006 0.008 ± 0.00421 

0.1387 0.1385 

0.55 fm 0.663 ± 0.006 22 

r^f^ 0.25 (dimensionless) 0.310 ± 0.003 



Table 3: Comparison between the calculated values of low-energy pion observables and ex- 
periment or, in the case of (— (^9)iGev)^^^ i^u/di the values estimated using other the- 
oretical tools. Each set in Eq. (69) yields the same calculated values. Aqcd = 0.234 GeV. 



magnitude of the pion form factor at large-g^. The two sets in Eq. (37) delimit reasonable 
boundaries and illustrate the model dependence in the result. In the chiral limit: /° ~ 
0.090 GeV. 

The quark propagator obtained with these parameter values is pointwise little different 
to that obtained in Ref . [23] . One gauge of this is the value of the Euchdean constituent 
quark mass. Here M^/^ = 0.32 GeV whereas M^f/^ = 0.30 GeV in Ref. [23]. It is also 
qualitatively similar to the numerical solution described in Sec. 4.1)^ where M,^^ = 
0.56 GeV. Indeed, the results are not sensitive to details of the fitting function: fitting 
with different confining, algebraic forms yields S{p) that is pointwise little changed, and 
the same results for observables. The earher parametrisation has been apphed more 
widely, as reviewed in Ref. [24] , and Table 4 summarises the results. 

In the calculations fT^r.,^ is 20% too small. This discrepancy cannot be reduced in impulse 
approximation because the nonanalytic contributions to the dressed-quark-photon vertex 
associated with tt-tt rescattering and the tail of the p-meson resonance are ignored. It 
can only be eliminated if these contributions are included. This identifies a constraint on 
realistic, impulse approximation calculations: they should not reproduce the experimental 
value of Jt^Tt^ to better than ~ 20%, otherwise the model employed has unphysical degrees- 
of-freedom. 

The pion form factor calculated using Eqs. (56) and (57) with (69) is compared with 
available data in Figs. 15 and 16. It is also compared with the result obtained in Ref. [23], 
wherein the calculation is identical except that the pseudovector components of the pion 
were neglected. Figure 15 shows a small, systematic discrepancy between the calculations 
and the data at low q^^ which is due to the underestimate of r^^ in impulse approximation.''" 
The results obtained with or without the pseudovector components of the pion Bethe- 
Salpeter amplitude are quantitatively the same, which indicates that the pseudoscalar 
component, is dominant in this domain. 

The increasing uncertainty in the experimental data at intermediate is apparent in 
Fig. 16, as is the difference between the results calculated with or without the pseu- 

tjust as in the present calculation, f-^TT^ = 0.25 in Ref. 23. However, the mass-scale is fixed so that 
/tt = 0.084, which is why this result appears to agree better with the data at small-g^: r^r is larger. 
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Calculated "Experiment" 





0.0924 GeV 


0.0924 ± 0.001 


fK 


0.113 


0.113 ± 0.001 




0.139 


0.138 


rriK 


0.494 


0.494 


1 11^ -17-2 
1 GcV 


0.0045 


0.008 


1 CjCV 


0.113 


0.1 ~ 0.3 


m/lGeV2 


n 9zL7 


u.zou m u.uuo 




u.oo im 


U.DDO ± U.UUD 




n /IQ 25 


n c;8Q 4-0 0/11 




-U.UZU IIli 


-U.U04 It U.U/u 




0.50 (dimensionless) 


n c;n/i -i- n ni q 
U.0U4 It u.uiy 




1.U4 


1 (Anomaly) 




u.iy 


U.ZD it u.uo 






n 098 4- 01 9 
-U.UZo It u.u±z 


2aS - 5a?. 


65 


n -L n 19 




0.038 


0.038 ± 0.002 


a° 


0.0017 


0.0017 ± 0.0003 


a| 


-0.00029 






1.22 


1.22 ± 0.01 




0.87 


0.88 ± 0.06 



Table 4: Summary of results obtained using the parametrisation of S{p) introduced in 
Ref. [23], which differs little from that specified by Eqs. (56) and (57). aj are tt-tt scat- 
tering lengths, whose experimental values are discussed in Ref. [28], and F^'^(4m^) is the 
value of the 771 tttt transition form factor at the softest point kinematically accessi- 
ble. The citations for the calculated results specify the article in which the annotated 
study is described. The "experimental" values of the current-quark masses and the quark 
condensate are estimates made using other theoretical tools: see Table 3. 



dovector components of the pion Bethe-Salpeter amplitude. These components provide 
the dominant contribution to ^^(q'^) at large pion energy because of the multiplicative 
factors: 7 • P and ^ ■ kk ■ P, which contribute an additional power of in the numera- 
tor of those terms involving F^, FG and relative to those proportional to E. Using 
the method of Ref. [23] and the model-independent asymptotic behaviour indicated by 
Eq. (63) one finds 

Fn[q ) oc — ; (70) 

y J -IT 

i.e., q^FT,{q^) ^ const., up to calculable Ing^-corrections. If the pseudovector components 
of are neglected, the additional numerator factor of is missing and one obtains 
q^FT^{q^) fa const. 

With this model the behaviour identified in Eq. (70) becomes apparent at > 2 M^y. 
This is the domain on which the asymptotic behaviour indicated by Eq. (63) is manifest. 
The calculated results, obtained with the two sets of parameters in Eq. (37), illustrate 
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Figure 15: Calculated pion form factor compared with data at small q^. The data are 
from Refs. [22] (crosses) and [30] (circles). The sohd hne is the result obtained when the 
pseudovector components of the pion Bcthc-Salpctcr amphtude are included, the dashed- 
hne when they are neglected. On the scale of this figure, both parameter sets in Eq. (69) 
yield the same calculated result. 



the model dependent uncertainty: 

q^FJq^) ~ 0.12 -0.19 GeV^ (71) 

^ ' q2^10-15GeV2 ^ ' 

It arises primarily because the model allows for a change in one parameter to be com- 
pensated by a change in another. In this example: 6f > 63 but &o + ^2 — ^0 + ^2 i 
and 6^63 = 6f 6f, which allows an equally good fit to low-energy properties but alters 
the intermcdiatc-g^ behaviour of F.,^{q^). As emphasised, in a solution of Eq. (7) these 
coefficients of the 1/p^ and 1/p® terms are correlated and such compensations cannot 
occur. 

As a comparison, evaluating the leading-order perturbative-QCD result with the asymp- 
totic quark distribution amplitude: 4)as{x) :— i/3/7rx(l — x), yields 

q'F^iq') = 16 7r/Xg2) ^ 0.15 GeV^ (72) 

assuming a value of a(g^ ~ 10 GcV^) ~ 0.3. However, the perturbative analysis neglects 
the anomalous dimension accompanying condensate formation; e.g., Eqs. (19)- (22) are 
not satisfied in Ref. [32]. 

5.2) Electroproduction of Vector Mesons. 
There is an extensive body of literature describing Pomeron phenomenology, all derived 
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Figure 16: Calculated pion form factor compared with the largest data available: 
diamonds - Ref. [30]; and circles - Ref. [31]. The sohd lines are the results obtained when 
the pseudovector components of the pion Bethe-Salpeter amplitude are included (lower 
line - set A in Eq. (69); upper line - set B), the dashed-line when they are neglected. 

from the observation that the total cross section in high-energy scattering: p-p, p-p, vr^-p, 
7-p, etc., is forward-peaked and rises slowly with ^Js. This is illustrated^^ in Fig. 17 and 
can be described by a Pomeron-exchange model of the quark-nucleon interaction with 
the following features: 

1. The quark-Pomeron coupling is q^{p2)^li(l^{Pi), where := /3/ 7^ with /3/ a flavour- 
dependent coupling constant. It is the only flavour-dependence in the interaction. 

2. The Pomeron "propagator" is characterised by a Regge trajectory: 

^(s,^) := (ais)"°+"i* (73) 

with ao > 0, which ensures the increase with s, and the Pomeron- nucleon coupling 
is 3/?„Fi(i), where Fi{t) is the Dirac form factor of the proton. 

3. The interaction is used in impulse approximation so that, for example, the ttN — > 
ttN interaction is completely described by 

{P;p2m'^\T^N-.7TN\q;Pims) := A^(g, P) SPu/dFi{t) G{s, t) Um'^ip2)l tiUm.ipi) , (74) 

where itms(Pi) is a nucleon spinor and A^(gf, P) :— 2AjJ(g, P) + 2A^(g, P) with 

Al{q,P) := N.tiD J ^5„(A;_+)F,(A;o-)S'rf(A;__)f,(A;)5„(A;++)/?/i7^ (75) 
describing the interaction of the Pomeron with the /-quark in the pion. 
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Figure 17: Total and elastic cross sections for p-p and p-p scattering. The slow increase 
of the total cross section with a/s at high energy is obvious. 

The parameters: Oq, ai, in this model were fixed by requiring a good description of 
7iN and KN scattering, and this is achieved with 

Pu/d = 2.35 GeV-2, = 1.50 GeV-^, = 0.10, ai = 0.33 GeV-^ . (76) 

In the diffractive regime the process e~ p —>■ e^' p'V, where V = p, (p, tp, is also expected 
to proceed via soft-Pomeron exchange and the model introduced above can be applied 
directly. The matrix element is 

{p2m2;kXp\Jf,\pimi) = 2f3ftf,^x{q,k) el''{k)Gv{w'^,t)3(3uFi{t)u{p2)^uu{pi) , (77) 

depicted in Fig. 18, where the 7P p transition form factor is 
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Figure 18: />meson electroproduction matrix element. oJ^ :— —{q — P/2 :— 

-{P+P2?- 



t,uxiq, P) = 3eo / ^ tr \s{k + \p) x (78) 
ri{k + \p,k-q + \p)S{k-q + \P) 7. S{k -\P)V^{k- -P) } . 

The unknown quantity in Eq. (78) is the vector meson Bethe-Salpeter amphtude, Vi,{k] —P) 
In the absence of a solution of the associated Bethe-Salpeter equation, an oft used and 
phenomenologically efficacious procedure is to parametrise the amplitude in a manner 
similar to that employed for the pion in Sec. 5.1): 

where Ny is fixed via the canonical normalisation condition: clearly, P ■V{k] P) = 0. The 

parameters are 

a(GeV) 6(GeV) c 
p 0.40 0.008 125.0 , . 

0.45 0.6 0.3 ^ ^ 

ijj 1.10 0.0 0.0 

which were fixed ^° by requiring the fit to the dimensionless coupling constants in Eq. (81). 

dp^e+e- dp^TT+n- 9cf>^e+e- 9<j)^KK 9ip^e+e- 

Theory 4.6 6.8 12.7 3.9 11.5 (81) 

Experiment 5.0 6.1 12.9 4.6 11.5 

The experimental values of the decay constants were calculated from the widths in 
Ref. [21], and the fit is acceptable given the simphcity of the Ansatz for V^, which in- 
cludes only one of the eight scalar functions necessary to completely specify a vector 
meson bound state. At this point there are no free parameters in the calculation of the 
electroproduction cross sections. 



32 




q' (GeV") 



Figure 19: p-meson electroproduction cross section aiW = 15 GeV: solid line - calculated 
result; long-dashed line - result if fhu/d — iOfhu/d', dash-dot line - result if rhu/d 
25fhu/d — frig. The data are: circles - Ref. [33]; squares - Ref. [34]; triangles - Ref. [35]. 

Figure 19 depicts the Q^-dcpendence of the p-meson electroproduction cross section and 
the magnitude is a prediction. There is complete agreement on the entire range of acces- 
sible Q^, with the large behaviour: l/Q^, which becomes evident at ~ 1-2 GeV^. 
Below that point the nonperturbative character of the dressed-quark propagator domi- 
nates the evolution of the cross section. It is important to observe the prediction that, 
the larger the current-quark mass of the constituents, the larger the value of at which 
the asymptotic behaviour is manifest. 

The calculated ^-meson electroproduction cross section is depicted in Fig. 20. It is in 
excellent agreement with Refs. [36] and [37], which used a nucleon target, as opposed to 
Ref. [35], which averaged over variety of nuclear targets. As could be anticipated from 
Fig. 19, the onset of the asymptotic 1/Q^ behaviour is pushed to larger-Q^ for the 0-meson 
because the current-quark mass of the constituents, the s-quark, is larger. 

In calculating the ■^-meson electroproduction cross section a very simple form was used 
for the c-quark propagator: 

Sc{k) := \ i-ij -k + mc) T{1 + jml) (82) 

with vric = 1.2 GeV (~ mj'^'^^ in Eq. (10)). This and the simple form of the ■?/'-meson 
Bethe-Salpeter amplitude anticipate the successful application of DSEs to heavy-meson 
observables. ^ The calculated cross section is depicted in Fig. 21. The striking prediction, 
confirmed by recent data, is that although two-orders of magnitude smaller than the p- 
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Figure 20: 0-meson electroproduction cross section atW = 70 GeV: solid line; the dashed 
hne is the p-meson result for comparison. The data are: triangles - Ref. [35]; circles - 
Ref. [36]; squares - Ref. [37]. 

meson cross section at the photoproduction point, the ■i/'-meson cross section is equal to 
the 1^ = 100 GeV, p-meson cross section at = 15 GeV^. This is because the large 
c-quark mass shifts the onset of the asymptotic l/Q'^-behaviour to larger-Q^. 

6. Finite Temperature and Chemical Potential. 

As we have seen, at zero temperature and chemical potential the low-energy and small- 
behaviour of QCD is characterised by confinement and DCSB. The internal scale that 
marks the boundary between small and large energy in QCD is ~ Aqcd- As the energy 
and/or momentum transfer increases, QCD is characterised by asymptotic freedom; i.e., 
the coupling evolves 



and quarks and gluons behave as weakly interacting, massless particles in high-energy 
and/or large-Q^ processes. 

The study of QCD at finite temperature and baryon number density proceeds via the 
introduction of the intensive variables: temperature, T; and quark chemical potential, 
fi. These are additional mass-scales, with which the coupling can run and hence, for 
T ^ Aqcd and/or fi ^ Aqcdj tts(Q^ = 0,T, yu) ~ 0. It follows that at finite temperature 
and/or baryon number density there is a phase of QCD in which quarks and gluons are 
weakly interacting, irrespective of the momentum transfer ^'^; i.e., a quark-gluon plasma. 
Such a phase of matter existed approximately one microsecond after the big-bang. In 



Q;s(g',T = = //) 








(83) 
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Figure 21: i/j-meson electroproduction cross section aX W = 100 GeV: solid line; the 
dashed line is the p-meson result at 1^ = 15 GeV for comparison. The data are from 
Refs. [38,39]. 

this phase confinement and DCSB are absent and the nature of the strong interaction 
spectrum is qualitatively different. The contemporary expectation for the position of the 
phase boundary in the (/i, T)-planc is illustrated in Fig. 22. 

The path followed in the transition to the plasma is important because it determines some 
observational consequences of the plasma's existence. For example, '^^ the time-scale for 
the expansion of the early universe: ~ 10~^ s, is large compared with the natural time-scale 
in QCD: 1/Aqcd ~ 1 fm/c ~ 10~^^ s, hence thermal equilibrium is maintained throughout 
the QCD transition. Therefore, if the transition is second-order, the ratio B :=baryon- 
numbcr/cntropy, remains unchanged from that value attained at an earlier stage in the 
universe's evolution. However, a first-order transition would be accompanied by a large 
increase in entropy density and therefore a reduction in B after the transition. Hence 
the order of the QCD transition constrains the mechanism for baryon number generation 
in models describing the formation of the universe, since with a second-order transition 
this mechanism is only required to produce the presently observed value of B and need 
not allow for dilution. In the absence of quarks, QCD has a first-order deconfinement 
transition, and with three or four massless quarks a first-order chiral symmetry restoration 
transition is expected. A current, primary question is what happens in the realistic case 
of two light quark fiavours? 

Nonperturbative methods are necessary to study the phase transition, which is charac- 
terised by qualitative changes in order parameters such as the quark condensate. One 
widely used approach is the numerical simulation of finite temperature lattice-QCD, with 
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Figure 22: The anticipated quark-gluon phase boundary in the (yUTv, T)-plane. "HG" - 
hadron gas, "QGP" - quark-gluon plasma. The nucleon chemical potential: fiN '■= 3/i.; 
i.e., three-times the quark chemical potential. "SPS" and "AGS" mark the points in the 
plane that it is estimated these facilities explore. 



the first simulations in the early eighties and extensive efforts since then. Here I focus 
on the application of DSEs. This is a new usage and much remains to be learnt: these 
are exploratory studies. One goal is to develop DSE models of QCD at finite-T and /i 
(QCD^) that can be used to check the results of numerical simulations, and be constrained 
by them. These models can then be employed to extrapolate into that domain presently 
inaccessible to lattice studies, such as finite chemical potential and the effects of T and /x 
on bound state properties, the latter of which are expected to provide the signatures of 
quark-gluon plasma formation in relativistic heavy ion collisions. 

Before discussing details it is interesting to provide a human scale for the temperatures 
and densities involved. The natural scale in QCD is Aqcd ~ 200 MeV and temperatures 
of this order of magnitude will be necessary to create the plasma. Aqcd ~ 10^° x room- 
temperature! It represents a temperature on the astrophysical and cosmological scale. 
Nuclear matter density po ~ 3 x 10^^ g/cm^ — 0.16 N/fm^ and this is more-than lO^^x 
the density of solid lead! The density at the core of a neutron star is expected to be 
approximately 4 po and it is densities on this order that are anticipated to be required 
for plasma formation. 

The expectation of the existence of a new phase of matter, the quark-gluon plasma, 
has led to the construction of a Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory. Construction is due to be completed in 1999. It will use counter- 
circulating, colliding 100 A GeV ^^^Au beams to generate a total centre-of-mass energy of 
~ 40 TcV, in an effort to produce an equilibrated quark-gluon plasma. It is anticipated 
to approach the quark-gluon plasma via a low baryon density route. Contemporary, fixed 
target experiments at the Brookhaven- AGS and CERN-SpS explore a high baryon density 
environment at much lower centre-of-mass energies. These experiments are crucial in 
developing the expertise necessary for operating detectors under RHIC conditions but they 
are not expected to "discover" the plasma. The CERN-SpS experiments have nevertheless 
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Figure 23: Mass spectra for inclusive e^-e~ pairs in 450 GeV p-Be collisions showing 
the data and various contributions from hadron decays. The shaded region indicates the 
systematic error on the summed contributions. 

produced some interesting results. 

One example is the "NA45-CERES" experiment, which studies e^-e~ pair production 
in heavy ion collisions, pairs leave the interaction region without interacting strongly 
and hence they are a probe of the early stages of the interaction. In Fig. 23 I illustrate the 
dilepton spectrum obtained in high-energy p-Be collisions. It is well described by known 
hadron decays. The same is true of p-Au collisions. However, this is not the case in S-Au 
collisions, as illustrated in Fig. 24. There the known hadron decays describe the data 
only for < 300 MeV, which is the region dominated by 7r° Dalitz decays. At higher 
energies the shape of the spectrum is different and shows a strong enhancement in the 
dilepton yield. Integrating over the region 0.2 < niee < 1.5 GeV the enhancement factor 
is 

5.0 ±0.7 Stat. ±2.0 syst. (84) 

The enhancement persists in Pb-Au collisions. What explanation can be offered? 

One model calculation^^ shows that this enhancement can be explained by a medium- 
induced reduction of the p-meson's mass, another that it follows from an increase in 
the p-meson's width. A decrease in the p-meson's mass is consistent with the QCD sum 
rules analysis of Ref. [48] but inconsistent with that of Ref. [49], which employs a more 
complex phenomenological model for the in-medium spectral density used in matching the 
two sides of the sum rule. In Ref. [49] there is no shift in the p-meson mass but a significant 
increase in its width. The consistency between Refs. [47] and [49] is not surprising since, 
in contrast to Ref. [48], they both rely heavily on effective Lagrangians with elementary 
hadron degrees-of-freedom. These are possibilities that can be explored using DSEs, 
which focusing on dressed-quark and -gluon degrees of freedom is an approach most akin 
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Figure 24: Mass spectra for inclusive e^- pairs in 200 GeV S-Au collisions showing 
the data and various contributions from hadron decays. The shaded region indicates the 
systematic error on the summed contributions. "^^ 



to Ref. [48]. A first attempt, summarised in Sec. 7.6), predicts a 15% suppression of 
nip at nuclear matter density but employs a model that is inadequate to address Vp. 

6.1) Notes on Field Theory. 

Equilibrium statistical field theory can be understood by analogy with equilibrium 
statistical mechanics. For a particle moving in a potential V the density matrix is given 
by the path integral 



p{x,x';T := 1/(3) := / Ps(r) exp <^ - / dr 

Jx{0)=x I JO 



^mx{Ty — V{x{t)) 



^5) 



where the underlined term is just the Lagrangian. All of the thermodynamic information 
about this system can be obtained from the partition function 

Z{T) := dxp{x,x,T); (86) 

for example, the pressure P = Tin Z{T)/V and the baryon density = {l/3)dP/dp. 

The equilibrium thermodynamics of a quantum field theory is also completely specified by 
a partition function, or generating functional. In the particular case of a self-interacting 
scalar field this is given by the functional integral: 

Z[T] := J nj,,e[o,/?] ^0(5, r) exp dr J d' x C"" {x; (f))^ , (87) 
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where C^{x; (f)) is the Euchdean Lagrangian density describing the interaction of r), 
whose boundary conditions are periodic: 

0(5, T = 0) = T = /3) . (88) 

The boundary conditions for fermions are antiperiodic: 

^(:r,r = 0) = -^(£,r = /5). (89) 

This difference in boundary conditions is the reason for the difference between the Mat- 
subara frequencies of fermions and bosons and hence why fermions acquire a screening 
mass at finite temperature. 

It is immediately obvious that the 0(4) invariance of the Euchdean theory is lost: at 
finite temperature (and/or chemical potential) the theory exhibits only an 0(3) sym- 
metry corresponding to spatial rotations and translations. This is why the formalism, 
necessarily used in lattice simulations, is only applicable to equilibrium systems - there is 

no generator of translations in time. One also notes from Eq. (87) that as T ^ cxo one 
dimension disappears completely and hence the corresponding {d— l)-dimensional theory 
is a candidate to describe the infinite-temperature limit of a d-dimensional theory. 

The finite temperature, free fermion Lagrangian density is 

^) = ^i^, r) (7 • a + 74 9, + m) ^{x, t) . (90) 

Introducing a Fourier decomposition: 

d^p 
(2^ 

antiperiodicity entails that the fermion Matsubara frequencies are 

Un = (2n+ l)7rr, n e Z (92) 

and the free fermion action is 



i^ix, t)^T Yl J jA-s m c^n)e^^-^+^'^"^ , (91) 



°° r d^p 

'S'|['0,'0]FVee = 2" J Tl^i^iP^^n) {il -P + ili^n + m) ij{p,UJn) . (93) 



n=— oo 



(27r) 

Prom this one identifies the finite temperature, free fermion propagator 

Sip) = ^ / . (94) 

l'^ ■ P + l'^4^UJn + m 

Analogous arguments, using the periodic boundary conditions, lead to an identification 
of the free boson propagator 

where Vln = 2 7inT. Having obtained the free particle propagators one can proceed to 
define a perturbation theory. As an example, in massless 0*^ theory the one-loop correction 
to the (f) propagator is 

^ f d^p 1 
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The sum can be evaluated: 



^ XI 02 I M2 = 7rTFF\ T + ^"-independent piece , (97) 

^fr'oo ^l + W \P\ exp(|pl/T) - 1 



to yield the Bose-Einstein factor. This is a source of problems: for large temperatures 

1 T 



(98) 



expdpl/T) - 1 \p\ 
and one can thereby encounter additional infrared divergences. 

These are particularly serious in QCD. To illustrate this^^ consider an {£ + l)-gluon- 
loop diagram and focus on the n = mode, which obviously yields the dominant infrared 
behaviour. The estimate is made easier if one neglects the tensor structure and notes that: 
from the vertices there is a factor of g^^ p"^^; the loop-sum-integral gives (T /p2g[0T] d^pY~^^', 
and the propagators, (p^ -|- m^)~^^ where m is a possible, dynamically generated mass- 
scale. A little thought and calculation shows that the net order of a given diagram in 
perturbation theory is 



1,2 e = 3 i>4 



Clearly, if m = the diagrams are infrared divergent for i > 2. The divergences may 
cancel when all diagrams of a given order are summed but that is difficult to verify. 
Suppose instead that the mass-scale m ~ ^fT, as does the Debye mass in QED, then 
no problem arises: at each order above i — 3 the diagrams are suppressed by powers of 
the coupling and a self-consistent calculation of the mass-scale is straightforward. This 
underlies the successful application of the method of "hard thermal loops" . However, 
suppose that m ~ g^T, which is the next possibility. In this case every diagram above 
i — 2 contributes with the same strength: g^, which presents a serious impediment to the 
application of perturbation theory! 

The introduction of a quark chemical potential modifies Eq. (93): 

Sp['tp,'tp]Free-^T ^ / — — ^'(p, C^n) («7 ' P + ^74 - 74/^ + "^) ^n) (100) 

n=-oo '' l^'^J 

SO that even the free Dirac operator is not hermitian and hence the QCD action is nec- 
essarily complex. As such it does not specify a probability measure, which precludes the 
straightforward application of Monte- Carlo methods in the evaluation of the partition 
function. However, the application of DSEs remains straightforward. The propagators 
and vertices are complex, so twice as many functions are required to represent them but 
otherwise there are no complications. Thus they provide a nonperturbative means of 
exploring this domain, which is presently inaccessible in lattice simulations. 

6.2) Some Lattice Results. 

Since the early eighties, as one branch of the extensive application of lattice methods 
in many areas of QCD, Monte-Carlo simulations have been used to estimate the finite 
temperature QCD partition function. These studies have contributed considerably to 
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Figure 25: Energy density and pressure for 2- light-flavour QCD on lattices with four 
temporal lattice sites, from Ref. [58]. 

the current understanding of the nature of the quark-gluon plasma. Due to the persis- 
tent limitation of computational power many analyses have focused on the pure gauge 
sector, which exhibits a first-order deconfinement transition at a critical temperature of 
Tc ~ 270 MeV. ^'^ In studying the chiral transition this commonly used quenched ap- 
proximation is inadequate because the details depend sensitively on the number of active 
(light) flavours. It is therefore necessary to include the fcrmion determinant. 

That is even more important when 7^ because the Dirac operator is not hermitian 
and thus the fcrmion determinant acquires an explicit imaginary part, in addition to 
those terms associated with axial anomalies. The QCD action being complex entails that 
the study of finite density is significantly more difficult than that of finite temperature. 
Simulations that ignore the fermion determinant at /x 7^ encounter a forbidden region, 
which begins at /i = m^/2,^^ and since — > in the chiral limit this is a serious 
limitation, preventing a reliable study of chiral symmetry restoration. The phase of the 
fermion determinant is essential in eliminating this artefact. 

QCD with dynamical quarks is a contemporary focus and for two flavours of light quarks 
the theory appears to exhibit a second-order transition at a critical temperature Tc ^~ ~ 
150 MeV. This is illustrated in Fig. 25, which shows a rapid change in the energy density 
in a small region around 150 MeV. For three or more light quark flavours the continuum 
theory is expected to have a flrst order chiral symmetry restoration transition. 

The quark condensate is an order parameter for chiral symmetry breaking, with its nonzero 
value at T = responsible for the pion mass being proportional to the square-root of the 
light current-quark masses. Its behaviour near the critical temperature is depicted in 
Fig. 26, and the rapid, qualitative change with increasing T is easily apparent. Very 
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Figure 26: Chiral condensate calculated^'' using staggered fermions and normalised to its 
zero temperature value. Only the Nf = results are extrapolated to zero quark mass. 



important is that, independent of the number of light-quark flavours, the condensate is 
unchanged for T < 0.9 Tg. It suggests that hadron properties are insensitive to T until 
very near the phase boundary. 

The simulations with dynamical fermions are still preliminary, and many uncertainties 
remain. For example, a review of recent results obtained with larger lattices and lighter 
quarks reports a significant lattice-volume-dependence for the critical exponents of the 
two light-flavour chiral symmetry restoration transition: the transition may even be first 
order! This might be an artefact of introducing lighter dynamical quarks, which drive the 
simulations to stronger coupling. The order of the transition with three and four flavours 
also remains unclear. With these uncertainties it is apparent that the lattice study of the 
phase transition will require further, even more computer-intensive simulations. 



7. DSEs at Finite T and 

The contemporary application of DSEs at finite temperature and chemical potential is 
a straightforward extension of the T — — /j, studies. The direct approach is to develop 
a finite-T extension of the Ansatz for the drcsscd-gluon propagator. The quark DSE can 
then be solved and, having the dressed-quark and -gluon propagators, the response of 
bound states to increases in T and /x can be studied. As a nonperturbative approach 
that allows the simultaneous study of DCSB and confinement, the DSEs have a signifi- 
cant overlap with lattice simulations: each quantity that can be estimated using lattice 
simulations can also be calculated using the DSEs. This means they can be used to check 
the lattice simulations, and importantly, that lattice simulations can be used to constrain 
their model-dependent aspects. Once agreement is obtained on the common domain, the 
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DSEs can be used to explore phenomena presently inaccessible to lattice simulations. 
7.1) Quark DSE. 

The renormalised dressed-quark propagator at finite- (T,/x) has the general form 

= -i7 -pcTAip, tOk) - iliL0k(7c{p, oJk) + crsip, tbk) (102) 

where ujk '■— oJk + i/J'- The complex scalar functions: A{p,u;k), B{p,u;k) and C{p,u;k) 
satisfy: 

J^{p, Uk)* = J^{p, oj-k-i) , (103) 

— A, B,C, and although not exphcitly indicated they are functions only of and o)^. 
The DSE for the renormalised dressed-quark propagator is 

S'"^(p, u}k) = Z^i^ -p^ Z2 {^4 Cbk + "^bm) + S'(p; ibk) , (104) 

where and Z2 are renormalisation constants, mbm is the bare mass, and the regularised 
self energy is 

S'(p,a'fc) = if -^2^(^,0;^) -Fi74cDfcS^(p,a'fc) + S^(p,a'fc) , (105) 

with 

f'^ 4: 1 

^'AP^^k) = / :i9^ Di^^{p- q,u}k-^^i)-7tr[Vjr'yi^S{q,0Ji)T^{q,0Ji;p,0Jk)^ 

Jl,q o 4 

where Va '■— —{Zf/p'^)i^ ■ p, Vb '■= Zi, Vc '■= — (■Z'i/a)fc)i74, Zf^ and Zi are vertex 
renormahsation constants, and Ii^q '■= T (i^' '^^^^ ^^^^ ^ mnemonic to 

represent a translationally invariant regularisation of the integral with A the regularisation 
mass scale. 

In Eq. (106) the Landau-gauge, finite-(T, //) dressed-gluon propagator has the form 



g'D^^ip, n) = P^M ^) ^f{p, ^) + P'iP) ^g{p, ^) (107) 



where 



0; fi and/or v — 4, 

Pi.(P) [s,-?f: = 1,2.3 . « 

E4 
Q=l PaPo 



Pi.m) -Pj,(p);/.,^^l,...,4. (109) 



A "Debye-mass" for the gluon appears as a T-dependent contribution to Ap. 
In renormalising we require that 

S~^{P,ujq) ° 2.2 = n + + (110) 

p -tu;q— c, 
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where C is the renormahsation point and is the renormahsed current-quark mass. This 
entails that the renormahsation constants are: 

Z^iCr^') = l-S>-',a'o)|g%g=,. , (111) 
Z,{CA') = l-S^(p,^o)lS;%g=,. , (112) 
= ^2mbm+ S^(p,'^o)|[j2%2=^2 , (113) 
and yields the renormahsed self energies: 

Tip, cl-fc) =^r + ^'Ap^ ^k) - S^(p, '^o) lf^2+^2^^2 , (114) 
where !F — A, B , C] — ^ — ic and — ttir- 

In studying confinement one cannot assume that the analytic structure of a dressed prop- 
agator is the same as that of the free particle propagator: it must be determined dy- 
namically. Indeed, one knows that the := (p, cjfc)-dependence of A and C is qualita- 
tively important since it can conspire with that of B to eliminate free-particle poles in 
the dressed-quark propagator. In this case the propagator does not have a Lehmann 
representation so that, in general, the Matsubara sum cannot be evaluated analytically. 
More importantly, it either complicates or precludes a real-time formulation of the finite 
temperature theory, which makes the study of nonequilibrium thermodynamics a very 
challenging problem. In addition, as we will see, this pfe-dependence of A and C can be 
a crucial factor in determining the behaviour of bulk thermodynamic quantities such as 
the pressure and entropy; being responsible for these quantities reaching their respective 
Stefan-Boltzmann limits only for very large values of T and fi. It is therefore important 
in any DSE study to retain A{pk) and C{pk), and their dependence on p^. 

7.2) Phase Transitions and Order Parameters. 

Phase transitions are characterised by the behaviour of an order parameter, {X), 
the expectation value of an operator. In the ordered phase of a system: {X) ^ 0, 
whereas in the disordered phase (X) = 0. A phase transition is first-order if (X) — > 
discontinuously, whereas it is second-order if (X) continuously. For a second- 
order transition, the length-scale associated with correlations in the system diverges as 
(X) — > and one can define a set of critical exponents that characterise the behaviour of 
certain macroscopic properties at the transition point. For example, in a system that is 
ferromagnetic for temperatures less than some critical value, Tc, the magnetisation, M, is 
an order parameter and in the absence of an external magnetic field M oc (T^ — T)^ for 
T ~ T~, where /3 is the critical exponent. At T = Tc the behaviour of the magnetisation in 
the presence of an external field, h 0"*", defines another critical exponent, 5: M oc hP^/^\ 
In a system that can be described by mean field theory these critical exponents are 

/3^^F = 0.5, 5^^ = 3.0. (115) 

The problem is that it can be difficult to identify the order parameter relevant to the 
discussion of a phase transition. 

One order parameter for the chiral symmetry restoration transition is well known - it is 
the quark condensate, defined via the renormahsed dressed-quark propagator: ^ 

-{,-,), := K Z,(i, A) £ , (116) 
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for each masslcss quark flavour, where the subscript "0" denotes that the scalar functions: 
Aq, Bq, Co, arc obtained as sohitions of Eq. (104) in the chiral hmit. The functions have 
an imphcit ^-dependence. An equivalent order parameter is 

A':=ReSo(p = 0,a;o), (117) 

which was used in Refs. [60-62]. Thus the zeroth Matsubara mode determines the char- 
acter of the chiral phase transition, a conjecture explored in Ref. [63]. 

What is an order parameter for deconfinement? 

In Sec. 2.3) I observed that the analytic properties of Schwinger functions play an impor- 
tant role in confinement. For illustrative simplicity, set = 0, the generalisation to 7^ 
is not difficult, and consider 

1 2 



Abo(x,t = 0) := T ^ / dppsm{px)aBo{p,u}n) (118) 

n=-oo ^'^•^ -"^ 

rp 00 

EASo(^)- (119) 



For a free fermion of mass M, (Jbo{P: <^n) = M/ (o;^ + p'^ + M^) and 



A^(a;) = Me-^V'-"+^' : (120) 



the n — term dominates the sum. In this case the "mass-function" 

d 

dx 



M{x;T) := (-In \A%^{x)\) - Vn^T^ + . (121) 



The most important observation is that for a free particle M{x, T) has a fixed, real value, 
which identifies the mass-pole in the propagator. It also exhibits the fermion "screening 
mass" = ttT, which becomes important for T ~ M/n. In the context of dynamical mass 
generation: M ~ . Since 450 MeV one anticipates that finite-T effects will 

become important at T ~ 150 MeV (or finite /i effects at ~ 450 MeV). For a boson of 
mass Mft, M{x; T) — M^: there is no screening mass. 

How does A°(x) behave if the dressed-propagator does not have a Lehmann representa- 
tion? An example is 

which has complex conjugate poles. In this case 

A?,(x) = e"^^ cos[hx\ ; (123) 

i.e., the Schwinger function oscillates and the mass-function has singularities, which is an 
unambiguous signal for the absence of a Lehmann representation and hence confinement! 

An order parameter for confinement is now obvious. Denote the position of the first zero 
in A^^(x) by which is inversely proportional to the distance of the poles from the real 
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Figure 27: —E{t) := InA(t) in QED3. Here the analogue of the mass-function is E'{t) 
and the difference between the unconfined theory: "+", and the confining theories is 
unmistakable. 



axis. Define kq '■— l/^oS then kq oc b and deconfinement is observed if, for some T — Tc, 
Ko{Tc) — 0: at this point thermal fluctuations have overwhelmed the confinement scale- 
parameter and the poles have migrated to the the real-axis. This criterion generalises 
easily to the case of ;U 7^ and to situations in which the drcsscd-propagator has an 
essential singularity rather than complex conjugate poles. It is also valid for both light 
and heavy quarks. 

An analogue of this criterion, with 

1 poo 

A{t) := — / dp.e^P^'asip^ 0,p,) , (124) 

ZTT J-00 

has been used to very good effect in an analysis of QED3 at T = 0. QED3 is confining 
in quenched approximation but not when massless fermions are allowed to influence the 
propagation of the photon. In that case complete charge screening is possible. Confine- 
ment is recovered in the theory if the fermion in the photon vacuum polarisation loop is 
massive. This application is summarised in Fig. 27. 

7.3) Illustration at {T ^0,n^0). 

As a first example I summarise a study ^° that uses a one-parameter, model drcsscd- 
gluon propagator. This parameter, m^, is a mass-scale that marks the boundary between 
the perturbative and nonperturbative domains, and its value, rrit — 0.69 GeV, was fixed in 
T — studies. The extension of the model to finite-T involves no additional parameters 
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and is defined with: ^f{p-, ^) '■— T^{p, fno) and ^g{p-, ^) '■— T^iP, 0); 



Vij), Q; m) 



27r 2. .3. ^ 1 _ e[-(P^+^^+-^)/(4-?)] 
— m^dond p) + .02^ — S 



(125) 



where = (8/3) tt^T^ is the perturbatively evaluated "Debye-mass" The quark DSE 
was solved using the rainbow approximation 



(126) 



I have discussed this truncation in Sees. 2.4) and 4), and here only note that in T = 
studies it has proven to be reliable in Landau gauge; i.e., at this level an efficacious 
phenomenology with a more sophisticated vertex Ansatz only requires a small quantitative 
modification of the parameters that characterise the small-A;^ behaviour of the dressed- 
gluon propagator. Using this truncation, mutually consistent constraints are Zf- — 
and Zx — Z2. 

The quark DSE was solved numerically with = 1.1 McV, C = 9.47 GeV. The 7 = 
fitting of mt and mji ensured a best x^-fit to a range of pion observables, yielding 



A = 92.4 
(92.4 ±0.3) 
al = 0.16 



m^r = 139.5 
(138.3 ±0.5) 

0,2 = -0.041 



r^N^ = 0.24 g^o^^ 
(0.31 ±0.004) (0.50 ±0.02) 



0.45 



0.02^ 



0.0022 



(127) 



(0.21 ±0.02) (-0.040 ±0.003) (0.038 ± 0.002) (0.0017 ± 0.0003) 



with the experimental values listed in parentheses.* The finite-T study reproduces these 
results to within 6% at T = 5 MeV, using the finite-T generahsations of the formulae in 
Ref. [67]: 

mlN'^ = {miiqq)c).; (128) 



(m^(gg)^)^ := SN^ f (ctbo - L^cr^ ± pV^ ± ct|] ) , 

Jk,p ^ '- -'' 



which vanishes linearly with m^; the canonical normalisation constant is 



4 2 



A 



k,p ^ 



i-2 



ujlaca'c +p'^aA(T'A + <JBa'j^] (129) 
with cr^ = daB{p^,u!k)/dp^, etc.; and the pion decay constant is obtained from 

f-A 



(130) 



^The influence of the Debye-mass on finite-T observables is quafitatively unimportant, even in the 
vicinity of the chiral symmetry restoration transition. The ratio of the coefficients in the two terms 
in Eq. (125) is such that the long-range effects associated with 5okS^{p) are completely cancelled at 



short-distances; i.e., for |xp ^ 1. 



•'■In Sec. 5.1) I discussed why TttNtt « 0.25 in impulse approximation. The tt-tt scattering lengths fitted 
in Ref. [67] were taken from Ref. [68]. 
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a 



X 1.1 GeV 

0.16 GeV 

Nl (0.18 GeV)2 

/,iV^ (0.15 GeV)2 
(m^(gg)) (0.15GeV)4 
0.12 GeV 
0.12 GeV 



0.33 
0.30 
1.1 
0.93 
0.92 
-0.11 
0.36 



Table 5: Parameters characterising the behaviour of the hsted quantities, fitted to 
q; (1 - T/T;)^, near = 150 MeV. 

Equations (128)-(130) were derived under the assumption that = ij^Bq. Some of the 
hmitations of this assumption were discussed in Sees. 4.2) and 5.1), and they are consid- 
ered further in Ref. [7]. It is quantitatively unreliable near the transition temperature, 
however, the qualitative behaviour of N^^ and /^r is the same, see Table 5. Only after 
these studies were completed was it understood that A?^ provides the best approximation 
to the leptonic decay constant when = ^Ts-Bq is assumed. 

The calculated T-dependence of the chiral symmetry and deconfinement order parameters 
is depicted in Fig. 28. The curves in the figure, fitted on T G [120, 150] MeV, are of the 
form a {1 — T/T^)^ with T^. ^ 150 MeV and a, j3 given in Tabic 5. The transitions are 
coincident and second-order with (3x = Pkq, within errors: ~ 10%. This estimate of j3x 
is not a mean field value and it agrees with a lattice estimate: ®^ (3^'^^ = 0.30 ± 0.08. It 
has been argued that two-hght-flavour QCD is in the universality class of the N — A 
Heisenberg magnet, for which = 0.38 ± 0.01 and both the DSE and lattice results 
are broadly consistent with this value. However, neither of these estimates of f3 survives 
more exhaustive study, ^^''^^ and the most recent analyses ^'"^^ suggest that in DSE models 
whose long-range part is described by the regularised singularity in Eq. (125) the chiral 
symmetry restoration transition at finite-T is described by a mean- field value of p. 

The behaviour of pion observables calculated from Eqs.(128)-(130) is depicted in Fig. 29. 
/tt and are weakly sensitive to T for T < 0.7 T^, and this is also seen in lattice simu- 
lations; e.g., the quark condensate in Fig. 26 and /^r in Fig. 30. However, as T approaches 
, the mass eigenvalue in the pion Bcthc-Salpeter equation moves to increasingly larger 
values, as thermal fluctuations overwhelm attraction in the channel, until at T = 
there is no solution and /^r — 0. This means that the pion-pole contribution to the four- 
point, quark- antiquark correlation function disappears; i.e., there is no quark- antiquark 
pseudoscalar bound state for T > . That may have important consequences for a wide 
range of physical observables, if borne out by improved studies; e.g., such T-dcpcndcnce 
for /jr and m^r would lead to a 20% reduction in the tt — > iiU/^ decay widths at T a; 0.9 T^"^. 

7.4) Complementary study at {T — 0, /i ^ 0). 
The finite-// behaviour of the same model has also been explored. The dressed-gluon 
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Figure 28: The order parameters for chiral symmetry restoration (A'(T), diamonds) and 
deconfinement {k,o{T), circles) both vanish at Tc — 150 MeV. The parameters for the 
fitted curves are presented in Table. 5. 
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propagator has the simple form* 

9'D,,{k) = U,.-:^)^; (131) 




G(P) 1R r 1 _ p-[feV(4m?)] 



;i32) 



and the rainbow approximation is used again. Neither the dressed propagator nor vertex 
have explicit /x-dependence, which can arise through quark vacuum polarisation insertions. 
As such they may be inadequate at large values of /i, particularly near any critical chemical 
potential. However, in the absence of finite-// studies of these quantities, the exploration 
of such models is useful, and one can assess the results obtained in the light of existing 
experiments and related theoretical studies. 

The renormalised dressed-quark propagator is 

S{P[,,]) ■■= -il ■ P(^a{P[,,]) -^74^Mf^c(PM) + > (133) 

where p[^] := {p,uj[^]), with uj[^] := + ifi. The quark DSE and the renormalisation 
conditions are similar to those discussed in the previous section, and the equation has two 
qualitatively distinct solutions: a chirally symmetric Wigner-Weyl mode, characterised 
by Bq = 0; and a confining Nambu-Goldstone mode, characterised Bq ^ 0. 

To explore the possibility of a phase transition one calculates the relative stability of 
the different phases, which is measured by the difference in pressure between them. The 
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Figure 29: Temperature dependence of the pion mass (m^(T), circles) and pion weak- 
decay constant (/^(T), diamonds). 



pressure is obtained directly from the partition function, Z: it is the sum of all vacuum- 
to-vacuum transition amplitudes. In "stationary phase" approximation, the partition 
function is given by the tree-level auxiliary-field effective action and the pressure is: 



P[S] :- 



T 
V 



T r 

- |TrLn 



Tr [J:S] 



(134) 



It is a functional of S{p[^]). In the absence of interactions E = and Eq. (134) yields 
the free fermion partition function. [Additive gluon contributions cancel in the pressure 
difference and are neglected.] The contribution of hadrons and hadron-like correlations to 
the partition function are neglected in Eq. (134). At the level of approximation consistent 
with Eq. (134) these terms arc an additive contribution that can be estimated using 
the hadronisation techniques of Ref. [74]. After a proper normalisation of the partition 
function; i.e., subtraction of the vacuum contribution, they are the only contributions to 
the partition function in the confinement domain. They are easy to calculate and are 
considered no further here as they are not a significant influence on the position of the 
phase boundary. 



The pressure difference is 



(135) 

which defines a /x-dependent "bag constant". In Eq. (135), A and C represent the 
solution of Eq.(104) obtained when Bq = 0; i.e., when DCSB is absent. This solution 
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Figure 30: Temperature dependence of the pion weak-decay constant on a 32^ x 8 lattice. 



exists for all fj,. B{ii) is plotted in Fig. 31. It is positive when the Nambu-Goldstone phase 
is dynamically favoured; i.e., has the highest pressure, and becomes negative when the 
Wigner pressure becomes larger. The critical chemical potential is the zero of i.e., 
/ic = 375 MeV. This abrupt switch from the Nambu-Goldstone to the Wigner- Weyl phase 
signals a first order transition. 

The order parameter for chiral symmetry restoration is that given in Eq. (117), while the 
confinement order parameter at {T — 0, /i ^ 0) is derived from 

1 r°° 

A5(r) := — / dp.e'P'^-aBoip^ 0,^;^) , (136) 

an analogue of Eq. (118). For a free, massive fermion ^^(p = 0,c<;[^]) = M/(c<;^^j + M^). 
This function has poles at pi = —{M±fi)'^, which are associated with the /^-induced offset 
of the particle and antiparticle zero-point energies, and 

A5(r) = ^e-(^-'^)-^(M-/.), (137) 

which is positive-definite and monotonically decreasing. In contrast, as observed above, 
for a Schwinger function with complex-conjugate p^-poles, As{t) has zeros at r > 0. 

The /^-dependence of the order parameters for chiral symmetry restoration and deconfine- 
ment is depicted in Fig. 32. The chiral order parameter increases with increasing chemical 
potential up to /ic, with X{iic)/X{0) ~ 1.2, whereas k{ii) is insensitive to increasing /i. At 
/ic they both drop immediately and discontinuously to zero, as expected of a first-order 
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Figure 31: B{n) from (135); B{fi) > marks the domain of confinement and dynamical 
chiral symmetry breaking. The zero of B{fi) is fic = 375 MeV. B{0) = (0.104 GeV)'', 
which can be compared with the value ~ (0.145 GeV)^ commonly used in bag- like models 
of hadrons.^^ 

phase transition. The increase of the chiral order parameter with is a necessary conse- 
quence of the momentum dependence of the scalar piece of the quark self energy, B{p[^]), 
as is easily seen in Ref. [61] and in Sees. 7.5) and 7.6). The vacuum quark condensate 
behaves in qualitatively the same manner as X. 

The behaviour of 771^^ and is illustrated in Fig. 33. One observes that although the 
chiral order parameter increases with /i, decreases slowly as increases. This slow fall 
continues until /i 0.7 /ic, when m,r(A*)/"T.7r(0) ~ 0.94. At this point m.j^ begins to increase 
although, for /i < /ic, i^-nilA does not exceed m^(0). This precludes pion condensation, 
in qualitative agreement with Ref. [76]. The behaviour of results from mutually 
compensating increases in {ra'R{<iq)(^)-K and A^^. This is a manifestation of the manner in 
which dynamical chiral symmetry breaking protects pseudoscalar meson masses against 
rapid changes with fi. The pion leptonic decay constant is insensitive to the chemical 
potential until ji ^ 0.7 /ic, when it increases sharply so that fTv{lJ'c)/f-^il^ = 0) ~ 1.25. 
The relative inscnsitivity of rrij^ and f^^ to changes in fi, until very near mirrors the 
behaviour of these observables at finite-T. For example, it leads only to a 14% increase 
in the tt — > /xi/ decay width at ~ 0.9 Hc- The universal scaling conjecture of Ref. [77] is 
inconsistent with the anticorrelation we observe between the //-dependence of /tt and mjr. 

Comparing the /i-dependence of and with their T-dependence, one observes an 
anticorrelation; e.g., at /t = 0, /^r falls continuously to zero as T is increased towards 
Tc ~ 150 MeV. This too is a necessary consequence of the momentum-dependence of 
the quark self-energy. In calculating these observables the natural dimension is mass- 
squared, and their behaviour at finite T and /i is determined by Re(a;p^]) ~ [tt^T^ — /f^]. 
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Figure 32: The order parameters for chiral symmetry restoration [X, diamonds] and 
deconfinement [k, circles] . = 375 MeV. 



where the T-dependence arises from the introduction of the fermion Matsubara frequency: 
P4 — (2n + l)7rT. Hence when such a quantity decreases with T it will increase with /i, 
and vice- versa. This is elucidated in Sees. 7.5) and 7.6), and in Ref. [50]. 

The confined-quark vacuum consists of quark-antiquark pairs correlated in a scalar con- 
densate. Increasing increases the scalar density: {—{qq))- This result is an expected con- 
sequence of confinement, which entails that each additional quark must be locally paired 
with an antiquark thereby increasing the density of condensate pairs as fi is increased. 
For this reason, as long as /i < /j.^, there is no excess of particles over antiparticles in the 
vacuum and hence the baryon number density remains zero; i.e., p^'^ = 0, Vfx < {ic- 
This is just the statement that quark-antiquark pairs confined in the condensate do not 
contribute to the baryon number density. 

The quark pressure, P^+'^f//], can be calculated, see Sec. 7.5), and one finds that af- 
ter deconfinement it increases rapidly, as the condensate "breaks- up", and an excess of 
quarks over antiquarks develops. The baryon-number density, p^"^ — {l/3)dP'^~^^/dfj,, 
also increases rapidly, with 

p^+'^(/i^2^,) ~3po, (138) 

where po = 0.16 fm^'^ is the equilibrium density of nuclear matter. For comparison, the 
central core density expected in a 1.4 Mq neutron star is 3.6-4.1 po- Finally, at // ~ 5//c, 
the quark pressure saturates the ultrarelativistic limit: P"+'^ = p^/(27r^), and there is a 
simple relation between baryon-density and chemical-potential: 

Pb"^'^/^) = J ^ , Vp > 5pe , (139) 

3 TT^ 

SO that p^'^'^'' {hpc) ~ 350 Po- Thus the quark pressure in the deconfined domain over- 
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Figure 33: Chemical potential dependence of the pion mass [m^, circles] and pion leptonic 
decay constant [/^, diamonds]. 

whelms any finite, additive contribution of hadrons to the equation of state, which antici- 
pating this was neglected in Ref. [62]. This discussion suggests that a quark-gluon plasma 
may be present in the core of dense neutron stars. 



7.5) Simultaneous study of {T ^ 0, /i ^ 0). 

This is the most difficult problem and the most complete study to date employs a 
simple Ansatz for the dressed-gluon propagator: 

g'D,,{p, n,) = (^6,, - J^^^ 27r3 ^ 4o , (140) 

which exhibits the infrared enhancement suggested by Ref. [6]. As an infrared- dominant 
model that does not represent well the behaviour of Q^) away from IjTI^ + ^ 0, 

some model-dependent artefacts arise. However, there is significant merit in its simplicity 
and, since the artefacts are easily identified, the model remains useful as a means of 
elucidating many of the qualitative features of more sophisticated Ansatze. 

With this model, using the rainbow approximation, the QCD^ gap equation, or DSE for 
the dressed-quark propagator, is^ 

S~'^ip, Uk) = Sq^{p, Uk) + \ri^lvS{p, LOk)lv ■ (141) 

A simplicity inherent in Eq. (140) is now apparent: it allows the reduction of an integral 
equation to an algebraic equation, in whose solution many of the qualitative features 
of more sophisticated models are manifest, as will become clear. In terms of the scalar 
functions introduced in Eq. (101), Eq. (141) reads 

jyW = + mB^ + (apI-t]'^ -m'^) -m (2r]'^ + + Apl) B , (142) 
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(143) 



Of particular interest is the chiral hmit, m = 0. In this case Eq. (142) reduces to a 
quadratic equation for B^pk), which has two quahtatively distinct solutions. The "Nambu- 
Goldstone" solution, for which 

B{Pk) = <( V'/ --^fe' ^<pI)<\ (144) 

otherwise 

Re(pi) < ^ 

. (145) 
otherwise , 





describes a phase of this model in which: 1) chiral symmetry is dynamically broken, 
because one has a nonzero quark mass-function, B{pk), in the absence of a current-quark 
mass; and 2) the dressed-quarks are confined, because the propagator described by these 
functions does not have a Lehmann representation. The alternative "Wigner" solution, 
for which 



,2 



Bip,) ^ , CiPk) = 1(^1 + ^1 + ||- J , (146) 

describes a phase of the model in which chiral symmetry is not broken and the dressed- 
quarks are not confined. 

With these two "phases", characterised by qualitatively different, momentum-dependent 
modifications of the quark propagator, this model can be used to explore chiral symmetry 
restoration and deconfincmcnt, and elucidate aspects of the method in such studies. 

In this model the relative stability of the different phases is measured by a (T, /i)- 
dependent "bag constant" , 

B{T,pl) := P[5ng]-/'[^w], (147) 

where S'ng means Eq. (101) obtained from Eq. (144) and ^w, Eq. (101) obtained from 
Eq. (146). As above, B{T, n) > indicates the stability of the confined (Nambu- 
Goldstone) phase and hence the phase boundary is specified by that curve in the (T, /x)- 
plane for which 

B(T,/x) = 0. (148) 



In the chiral limit 



B{T, /.) = rj' 2N,Nf- ^ dyy' |Re {2pf) - Re j - In \pfC{pi) 



(149) 



with: T — T/r), iJ, — /ji/r); l^ax is the largest value of I for which ouf^^^ < ^ + A*^ and this 

'^Lax ~ Pi ~ (y^ + ^A*)- '^(^' A*) depicted in Fig. 34 and 



also specifies 
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Figure 34: B{T,fi) from Eq. (149); B{T,fi) > marks the confinement domain. The scale 
is set by i3(0, 0) = (0.102 r/)^ = (0.109 GeV)^ r] = 1.06 GeV.^^ 



the critical line in Fig. 35. The deconfinement and chiral symmetry restoration transitions 
are coincident. 

For /i = the transition is second order and the critical temperature is T° = 0.159 rj, which 
using the value of ?7 = 1.06 GeV obtained by fitting the vr and p masses corresponds to 
T° = 0.170 GeV. This is only 12% larger than the value obtained in Sec. 7.3),^° and the 
order of the transition is the same. However, in the present case the critical exponent 
is (3 = 0.5. For any fi ^ the transition is first-order, as revealed by close scrutiny of 
Fig. 34. For T = the critical chemical potential is /i° = 0.3 GeV, which is ~ 30% smaller 
than the result in Sec. 7.4). One notes from Fig. 35 that Pc{T) is insensitive to T until 
T 0.3 T°. The discontinuity in the order parameters vanishes as /i — >■ 0. 

In the deconfinement domain, illustrated clearly in Fig. 35, the quarks contribute an 
amount 



P[S- 



w 



rp OO 



dyy"^ I In 



- 1 + Re 



1 



(150) 



to the pressure, which must be renormalised to zero on the phase boundary. Just as for 
free fermions, this expression is formally divergent and one must isolate and define the 
active, temperature-dependent contribution. This is difficult because, in general, C{pi) 
is only known numerically and hence it is not possible to evaluate P[5'w] analytically. A 
method for the numerical evaluation of Eq. (150) was developed in Ref. [61]. 
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Figure 35: The phase boundary in the (T,/i)-plane obtained from (148) and (149). The 
"structure" in this curve, apparent for small-T, is an artefact of the inadequate represen- 
tation of the quark-quark interaction in the ultraviolet by Eq. (140). 



Consider the derivative of the integrand in Eq. (150): 



1=0 



1=0 



(151) 



{y - fxf + ujj {y + + ujt 



V c{piY df 



In 



the absence of interactions C{pi) = 1, the second term is zero and 



2 °° 

rp 5Z 

^ 1=0 



jy - P? 



_{y-py + uf ^ {y + pY + uf 
where in this case e(y) = y and 



(y + ^^ 



r,\2 



(152) 



m = In 



1 -|- exp 



T 



+ In 



1 -|- exp 



(153) 



Appropriately inserting Eq. (152) for the parenthesised term in Eq. (150), and neglecting 
T-independent terms one obtains, 

'T roo 

(154) 



J" POO 

P[So] = v'N^Nf- dyy^I(y) 



7 



(155) 



which is the massless free particle pressure. 
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Figure 36: K{y,p,), which describes the nonperturbative modification of the free particle 
dispersion law, for p, — 0, 0.3, 0.6. By assumption, it is independent of T. 

To proceed in the general case the assumption is made®^ that the nontrivial momentum 
dependence of C{pi), which is manifest in all DSE-models of QCD^, acts primarily to 
modify the usual massless, free particle dispersion law. One evaluates the sum on the 
right-hand-side of (151) numerically and uses the form on the right-hand-side of Eq. (152) 
to fit a modified, T-independent dispersion law, e{y,Ji) = y + n{y,fi), to the numerical 
results. The existence of a K,{y, pi) that provides a good x^-fit on the deconfinement domain 
is understood as an a posteriori justification of the assumption. In Ref. [61] the relative 
error between the fit and the numerical results is < 10% on the entire T-domain. 

The calculated form of /t(?/, p) is depicted in Fig. 36; it only depends weakly on p. The 
form indicates a persistence of nonperturbative effects into the domain of deconfinement, 
evident in the nontrivial momentum dependence of C{pi) and its slow evolution to the 
asymptotic value C{pi) — 1. The effect of this is to generate a mass-scale in the massless 
dispersion law: k(0, 0) ^ 0.6 ~ 2/2^. This mass-scale is unrelated to the chiral- symmetry 
order parameter, X in Eq. (117), and is a qualitatively new feature of the study. For 
p > 5/i° the explicit mass-scale introduced by the chemical potential overwhelms the 
dynamically generated scale. 

Using this result, Eq. (150) becomes 

P[5w] = v'N^Nf- dyy'X{e{y,p)), (156) 

TT^ Jo 

and the quark pressure in this DSE-model of QCD^ is 

P,(T,/x) = eiV) {P[5w] - P[Sw]\av} (157) 

where V is the domain marked "Deconfined" in Fig. 35, 0{V) is a step function, equal 
to one for (T, /x) G V, and P[>S'w]|ax> indicates the evaluation of this expression on the 
boundary of P, as defined by the intersection of a straight- line from the origin in the 
(T, /x)-plane to the argument- value. It is plotted in Fig. 37, which illustrates clearly that 
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Figure 37: The quark pressure, Pq{T,fi), normalised to the free, massless (or Ultra- 
Relativistic) result, Eq. (155). 

in this model the free particle (Stefan-Boltzmann) limit is reached at large values of T 
and /i. The approach to this limit is slow, however. For example, at T ~ 0.3 ~ 2T^, 
or /i ~ 1.0 ~ 3/i°, Eq. (157) is only 50% of the free particle pressure, Eq. (155). A 
qualitatively similar result is observed in numerical simulations of lattice-QCD actions at 
finite-T. This feature results from the slow approach to zero with y of K{y, /i), illustrated 
in Fig. 36, and emphasises the persistence of the momentum dependent modifications of 
the quark propagator. 

With the definition and calculation of the pressure, Pq{T,n), all the remaining bulk ther- 
modynamic quantities that characterise the model can be calculated. As an example 
the "interaction measure": A := e — 3P, where e is the energy density, is plotted in 
Fig. 38. It is zero for an ideal gas, hence the name: A measures the interaction-induced 
deviation from ideal gas behaviour. This figure provides a very clear indication of the 
persistence of nonperturbative effects into the deconfinement domain, with a /i = max- 
imum of A ^ 0.2 P[So\ at T ^ 2T^ and a T = maximum of A ^ 0.3 P[So] at /i ^ 3/1^. 
Both Figs. 37 and 38 indicate that there is a "mirroring" of finite-T behaviour in the 
/i-dependence of the bulk thermodynamic quantities. 

7.6) 7c and p properties. 

The model discussed in the last section has been used ^° to study the (T, /i)-dependence 
of 71 and p properties, and to elucidate other features of the models described above 
that employ a more sophisticated Ansatz for the dressed-gluon propagator. In these 
applications its simplicity is particularly helpful. 
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Figure 38: The "interaction measure", A(T, /i), normalised to the free, massless result for 
the pressure, Eq. (155). 



To begin, consider the vacuum quark condensate, which in this model is 

-{QQ)=v'^fY. I dyy'Rel^^\-y^-ufy. (158) 

for T = = yU, {—{qq)) = ?7'^/(80 7r^) = (0.11 r/)^. In Fig. 39 one observes that {—{qq)) 
decreases with T but increases with increasing n, up to a critical value of fic{T) when 
it drops discontinuously to zero. These results are in qualitative and semiquantitative 
agreement with the (T = 0, fi ^ 0) and (T 7^ 0,/i = 0) studies described in Sees. 7.3) 
and 7.4). The increase with /i is also qualitatively identical to that observed in a ran- 
dom matrix theory with the global symmetries of the QCD partition function. {—(qq)) 
must increase with /i in the confinement domain because confinement entails that each 
additional quark must be locally-paired with an antiquark, thereby increasing the den- 
sity of condensate pairs. This vacuum rearrangement is manifest in the behaviour of the 
necessarily-momentum-dependent scalar part of the quark self energy, B{pk). 

In this model Eqs. (128)-(130) yield very simple expressions in the chiral limit; for 
example,'^ 

= v'^TY.\il + ^n'-^u.l-\4). (159) 

Characteristic in Eq. (159) is the combination /i^ — ujf, which entails that, whatever 
change f^, undergoes as T is increased, the opposite occurs as is increased. Without 

TThis is the expression for N"^ from Eq. (129), which provides provides a better approximation to the 
pion leptonic decay constant than Eq. (130) when one assumes Tt^{p] P) — 175 Bo{p^). 
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Figure 39: The quark condensate, Eq. (158), as a function of for a range of values of 
T. In all existing studies, in which the quark mass function has a realistic momentum 
dependence, it increases with n and decreases with T. At the critical chemical potential, 
ndT), {—{qq)) drops discontinuously to zero, as expected of a first-order transition. For 
/i = it falls continuously to zero, exhibiting a second-order transition at Tc(/x = 0) = 
0.16r/. 

calculation, Eq. (159) indicates that will decrease with T and increase with /i. This 
provides a simple elucidation of the results described above. Figure 40 illustrates this 
behaviour for m 7^ 0. The (T, yu)-dependence of rriT^, from Eq. (128), is also depicted 
in Fig. 40. It is insensitive to changes in fi and only increases slowly with T, until T 
is very near the critical temperature. As in Sec. 7.4), this insensitivity is the result of 
mutually cancelling increases in {mqq)^^ and /^r, and is a feature of studies that preserve 
the momentum-dependence of the confined, dressed-quark degrees of freedom in bound 
states. 

With 77 = 1.37 GeV and m = 30MeV, one obtains = 92MeV and = 140 MeV at 
T — — fj,. That large values of rj and m are required is a quantitative consequence of the 
inadequacy of Eq. (140) in the ultraviolet: the large-p^ behaviour of the scalar part of the 
dressed-quark self-energy is incorrect. This defect is remedied easily^ without qualitative 
changes to the results presented here. 

p-meson properties are more difficult to study: one must solve the vector-meson Bethe- 
Salpeter equation directly. As described above, the ladder truncation of the kernel in the 
inhomogeneous axial- vector vertex equation and the rainbow truncation of the quark DSE 
form an AV-WTI identity preserving pair. It follows that the ladder BSE is accurate for 
fiavour-nonsinglet pscudoscalar and vector bound states of equal-mass quarks because 
of a cancellation in these channels between diagrams of higher order in the systematic 
expansion illustrated in Fig. 7. 
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Figure 40: The pion mass, Eq. (128), and weak decay constant, Eq. (129), as a function 
of II for a range of values of T. m^r falls slowly and uniformly with /i [mT^{T — 0, fic) — 
0.95 m^(T = 0,ij, — 0)] but increases with T. Such a decrease is imperceptible if the 
ordinate has the range in Fig. 41. /^r increases with /i and decreases with T [fn{T — 
0, lie) = 1.51 MT = 0,11 = 0)]. 

A ladder BSE using the T = limit of Eq. (140) was introduced in Ref. [12]. It has one 
notable pathology: the bound state mass is determined only upon the additional specifi- 
cation that the constituents have zero relative momentum. That specification leads to a 
conflict with Eqs. (19)-(22), which follow from the AV-WTI, and is an artefact of imple- 
menting the delta-function limit discontinuously; i.e., these identities are manifest for any 
finite-width representation of the delta-function, as this width is reduced continuously to 
zero. In other respects this ladder BSE provides a useful quahtative and semi-quantitative 
tool for analysing features of the pseudoscalar and vector meson masses. For example, 
Goldstone's theorem is manifest, in that the vr is massless in the chiral limit, and also 
rises linearly with the current-quark mass. Further, there is a naturally large splitting 
between and rUp, which decreases slowly with the current-quark mass. 

To illustrate this and determine the response of rUp to increasing T and /i, the BSE of 
Ref. [12] was generalised to finite-(T, /x) as 

TMiPk; Pe) = - J Re {tm S{P, + ip,) TM{p^; Pi) S{p, - ip,) , (160) 

where P^ := (P, He) ■ The bound state mass is obtained by considering Pe=o and, in ladder 
truncation, the p- and a;-mesons are degenerate. 

The TT equation admits the solution 

r^{Po)^i5{iOi + i-pe2) (161) 
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Figure 41: Mp+ and function of fx for T = 0, 0.1. On the scale of this figure, 771.,^ 

is insensitive to this variation of T. The current-quark mass is m = 0.011 77, which for 
77 = 1.06 GeV yields Mp+ = 770 MeV and = 140 MeV at T = = 



and yields the mass plotted in Fig. 41. The mass behaves in qualitatively the same manner 
as in Fig. 40, from Eq. (128), as required if Eq. (160) is to provide a reliable guide. 
In particular, it vanishes in the chiral limit. 

For the p-meson there are two components: one longitudinal and one transverse to P. 
The solution of the BSE has the form 



74^, 



(162) 



where Opj^ labels the longitudinal and 9p_ the transverse solution. The eigenvalue equation 
obtained from Eq. (160) for the bound state mass, Mp±, is 



Re < 0-5 (a;. 



2 

0+ 



0+ 



4 p± 



X)1 = i- 



(163) 



The equation for the transverse component is obtained with [— ij-'o+ ~ (163). 
Using the chiral- limit solutions, Eq. (144), one obtains immediately that 

= i 77^, independent of T and (164) 

This is the T = = fi result of Rcf. [12]. Even for nonzero current-quark mass, Mp_ 
changes by less than 1% as T and /x are increased from zero toward their critical values. 
Its insensitivity is consistent with the absence of a constant mass-shift in the transverse 
polarisation tensor for a gauge-boson. 
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For the longitudinal component one obtains in the chiral limit: 

= \rf - 4(/.2 - ^T^) . (165) 

The characteristic combination [fj^ — n'^T^] again indicates the anticorrelation between 

the response of Mp^ to T and its response to yU, and, like a gauge-boson Debye mass, that 
M^_^ rises linearly with for /i = 0. The m 7^ solution of Eq. (163) for the longitudinal 
component is plotted in Fig. 41. As signalled by Eq. (165), Mp^ increases with increasing 
T and decreases as increases.'^ 

1 stated that contributions from skeleton diagrams not included in the ladder trunca- 
tion of the vector meson BSE do not alter the calculated mass significantly because of 
cancellations between these higher order terms. ^ This is illustrated explicitly in two cal- 
culations: Ref . [79] , which shows that the p — > nn — > p contribution to the real part of the 
p self-energy; i.e., the n-n induced mass-shift, is only —3%; and Ref. [80], which shows, 
for example, that the contribution to the a;-meson mass of the uj — > 37r-loop is negligible. 
Therefore, ignoring such contributions does not introduce uncertainty into estimates of 
the vector meson mass based on Eq. (160). 

Equation (163) can also be applied to the 0-meson. The transverse component is insensi- 
tive to T and /i, and the behaviour of the longitudinal mass, M<^+, is qualitatively the same 
as that of the p-meson: it increases with T and decreases with /j,. Using 77 = 1.06 GeV, 
the model yields M^± = 1.02 GeV for = 180 MeV at T = = /x. 

In a 2-flavour, free-quark gas at T = the baryon number density is ps — 2p^/{3n^) , by 
which gauge nuclear matter density, po = 0.16 fm~^, corresponds to p, — :— 260 MeV = 
0.245 7]. At this chemical potential the algebraic model yields 

Mp^po) ^ 0.75Mp4p^0), (166) 
M^+{po) ^ 0.85M^+(// = 0) . (167) 

The study summarised in Sec. 7.4),^^ indicates that a better representation of the ultra- 
violet behaviour of the dressed-gluon propagator expands the horizontal scale in Fig. 41, 
with the critical chemical potential increased by 25%. This suggests that a more realistic 
estimate is obtained by evaluating the mass at p'^ — 0.20 rj, which yields 

Mp+(p',) ^ 0.85M,+(/i = 0) , (168) 
M^+{p'q) ^ 0.90M^+(/x = 0) ; (169) 

a small, quantitative modification. The difference between Eqs. (166) and (168), and that 
between Eqs. (167) and (169), is a measure of the theoretical uncertainty in the estimates 
in each case. This reduction in the vector meson masses is quantitatively consistent with 
that calculated in Ref. [48] and conjectured in Ref. [81]. At the critical chemical potential 
for T = 0, Mp+ fa 0.65 Mp+{p = 0) and ^ 0.80 M^+{p = 0). 

This simple model of QCD^ preserves the momentum-dependence of gluon and quark 
dressing, which is an important qualitative feature of more sophisticated studies. Its 

t There is a 25% difference between the value of 77 required to obtain the T = = /i vahies of m,-,^ and 
fjr, from Eq. (128) and Eq. (129), and that required to give Mp± = 0.77 GeV. This is a measure of the 
quantitative accuracy of this algebraic model. 
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simplicity means that many of the consequences of that dressing can be demonstrated 
algebraically. For example, it elucidates the origin of an anticorrelation, found for a range 
of quantities, between their response to increasing T and that to increasing /x. 

Both {—{qq)) and decrease with T and increase with /i, and this ensures that 171^^ is 
insensitive to increasing /j, and/or T until very near the edge of the domain of confinement 
and DCSB. The mass of the transverse component of the vector meson is insensitive to 
T and ^ while the mass of the longitudinal component increases with increasing T but 
decreases with increasing This behaviour is opposite to that observed for {—{qq)) and 
and hence the scaling law conjectured in Ref. [81] is inconsistent with this calculation, 
as it is with others of this type. 

This study has two primary limitations. First, the width of the vector mesons cannot be 
calculated because the solution of Eq. (160) does not provide a realistic Bethe-Salpeter 
amplitude. Second, the calculation of meson-photon observables at T — = /j, only 
became possible with the determination ^ of the form of the dressed-quark-photon vertex. 
Its generalisation to nonzero- (T, /i) is a necessary precursor to the study of these processes. 

8. Closing Remarks. 

These lecture notes illustrate the contemporary application of Dyson-Schwinger equa- 
tions to the analysis of observable strong interaction phenomena, highlighting the positive 
aspects and successes. Many recent, interesting studies have been neglected; a calcula- 
tion of the electric dipole moment of the p-meson ®^ and an exploration of rj-r]' mixing 
among them. However, a simple enquiry of "http://xxx.lanl.gov/find/hep-ph" with the 
keywords: "Dyson-Schwinger" or "Schwinger-Dyson" , will provide a guide to other cur- 
rent research. 

In all phenomenological applications, modelling is involved, in particular, of the behaviour 
of the dressed Schwinger functions in the infrared. [The ultraviolet behaviour is fixed be- 
cause of the connection with perturbation theory.] This is tied to the need to make 
truncations in order to define a tractable problem. Questions will always be asked regard- 
ing the fidelity of the modelling. The answers can only come slowly as, for example, more 
is learnt about the constraints that Ward Identities and Slavnov- Taylor identities in the 
theory can provide. That approach has been particularly fruitful in QED,^ and already 
in the development of a systematic truncation procedure for the kernel of the quark DSE 
and meson BSE. '^'^ In the meantime, and as is common, phenomenological applications 
provide a key to understanding which elements of the approach need improvement: one 
must push and prod to find the weak links. 
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